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Abstract—This paper forms an introduction to a series dealing with the characterization of 
ion-exchangers both as regards their equilibria with electrolyte solutions and their kinetic 
behaviour. An outline is given of the development of ion-exchangers with reference to their 
principal field of application, the purification of water 


For the quantitative description of the behaviour of a column use will be made of an equation 
giving the differential material balance, and of data on equilibria and kinetics. When commercial 
products are examined it is, in general, found that the equilibria in which H*-ions play a part 
cannot be described by the known formulae, e.g., that of the law of mass action, while the deter- 
mination of their equilibria by known methods also meets with difficulties. 


The ion-exchange phenomenon is studied in the light of the theory of adsorptive percolation 
or chromatography. If there is an instantaneous establishment of equilibrium, the differential 
material balance may be used for demonstration how a chromatogram can alter its shape depending 
on the form of the isotherm. Mutatis mutandis, the same formulae as used for chromatography 
may be employed for ion-exchange. 


Experiments and applications, which will be discussed later, relate to the use of ion-exchangers 
in water purification practice 


Résumé— La présente documentation est une introduction a une série relative aux caractéristiques 
d'échangeurs d’ions tant en ce qui concerne leurs équilibres avec des solutions électrolytiques 
que leur comportement cinétique. On donne une vue générale du développement des échangeurs 


d'ions par rapport & leur domaine principal d’application, & savoir lépuration d'eau. 


Pour la description quantitative du comportement d'une colonne on utilise une équation 
représentant le bilan-matiére differentiel et des données relatives aux équilibres et a la cinétique. 
En examinant des produits de commerce on constate, en général, que léquilibre dans lequel 
les ions H* jouent un réle ne peut-@tre décrit moyennant les formules connues, par exemple, 
celle de la loi sur action de la masse, parce que la détermination de leur équilibre par les méthodes 
connues se heurte a des difficultés. 


On étudie le phénoméne de l'échange dions a la lumiére de la théorie de la filtration adsorptive 
ou la chromatographie. En cas d'un établissement instantané de léquilibre on peut utiliser 
un bilan-matiére differentiel pour démontrer comment un chromatogramme peut modifier sa 
forme en fonction de la forme de lisotherme. Les formules utilisées pour la chromatographie 
s’‘appliquent, mutatis mutandis, aussi & l'échange dions. 


Les expérimentations et les applications qui seront traitées plus tard portent sur l'emploi 
d’échangeurs d’ions pour l’épuration d'eau. 
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Zusammenfassung—Dieser Aufsatz stellt eine Einfiihrung zu einer Serié dar, die sich mit den 
charakteristischen Eigenschaften von Ionenaustauschern beschiiftigt und zwar sowohl hinsichtlich 
ihrer Gleichgewichte mit elektrolytischen Lésungen als auch ihres kinetischen Verhaltens. 
Ferner wird ein Uberblick gegeben tiber die Entwicklung von Ionenaustauschern in Bezug auf 
ihre hauptsichliche Anwendung, die Wasserreinigung. 


Fiir die quantitative Beschreibung des Verhaltens einer Kolonne wird von einer differen- 
ziellen Mengenbilanz und von den Werten des Gleichgewichts und der Kinetik Gebrauch gemacht. 
Bei der Untersuchung handel siiblicher Produkte stellt man fest, dass im allgemeinen die 
Gleichgewichte, in denen H*-Ionen eine Rolle spielen, nicht durch bekannte Gleichungen 
beschrieben werden kénnen, z.B. durch das Massenwirkungsgesetz, da auch die Bestimmung 
ihrer Gleichgewichte durch bekannte Methoden auf Schwierigkeiten stésst. 


Der Ionenaustausch wird unter dem Gesichtspunkt der adsorptiven Filtration oder Chromato- 
graphie untersucht. Sofern sich das Gleichgewicht augenblic klich einstellt, kann die differenzielle 
Mengenbilanz benutzt werden, um die Anderung der Form eines Chromatogramms in Abhingig- 
keit von der Form der Isothermen zu zeigen. Dieselben Gleichungen, die fir die Chromato- 
graphie benutzt werden, kénnen auch mutatis mutandis auf den Ionenaustausch angewendet 


werden. 


Experimente und Anwendugen, die spiter diskutiert werden sollen, beziehen sich auf den 
Gebrauch von Ionenaustauschern in der Praxis der Wasserreinigung. 


1 OBJECT OF THE INVESTIGATION 


JON-EXCHANGE, i.e. the process which 
equival nt amounts of ions are exchanged between 
a solid and a solution, was discovered about a 
century ago by THOMPSON [12] and Way (14). 
A few years later E1cunorn [5] established the 
reversibility of the process. Nearly 50 years 
after this discovery Gans [6] realized the first- 
technical application of ion-exchangers, water- 
softening, which is still their principal field. 

Up to about 20 years ago natural zeolites and 
synthetic silicates, for which Gans introduced the 
name permutites, were the only products used 
for this process. Gradually the permutites were 
superseded by ion-exchangers prepared from 
sulphonated coal [11] and LirBKNECHT 
[10]) and by completely synthetic products 
derived from phenosulphonic acid and formal- 
dehyde (Apam and Houmes [1]) which have the 
advantage of permitting regeneration with acids. 
The latter authors also developed the first 
anion-exchangers on the basis of polyamides and 
formaldehyde, while the range of ion-exchangers, 
as far as their functional groups are concerned, 


was completed by the strongly basic anion- 


exchangers (p’ALELIO [2]) which are capable 


even of binding weak acids such as silicic acid. 
As in other fields of technology, preparation 


and utilization of ion-exchangers has been far 
ahead of the theoretical knowledge. We may 
say that it was not until 1943, when de Vautt [13] 
advanced his theory of chromatography, that a 
basis was provided permitting a theoretical 
interpretation of the process. Of the problems 
attending the application of ion-exchangers, 
their characterization is one of the most important; 
unfortunately, this is still very incompletely 
understood. The cause of this is to be sought 
in the complicated character of the problem. 
For complete comprehension it is essential that 
information be available about the equilibrium 
and the kinetics of the process, and these data 
should be combined with the material balance 
of the column expressed in the two variables 
place and time. Since even in relatively simple 
cases the results are already extremely complicated 
from a mathematical point of view, it is obvious 
that kinetic investigations in a column are almost 
exclusively restricted to those processes of which 
the equilibria cither can be represented by a 
linear [4] or a Langmuir [9] isotherm or satisfy 
the law of mass action [3]. 

However, the results of this study of industrially 
applicable commercial products have shown 
that the above equilibria usually fail to fit the 
experimental results adequately. This also 
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appears from the literature dealing for example 
equilibrium measurements on 
exchangers in which H’*-ions involved 
(KrisHNaAMoorty and Overstreet [8]); it is 
just these equilibria which occur frequently in 


with cation- 


are 


connection with the use of ion-exchangers for 
water purification. 

Furthermore it will be seen that the customary 
procedures of determining the equilibria are not, 
types of 
exchangers. Therefore, another method has been 


in general, applicable to all ion- 
developed starting from the theory of chromato- 
graphy. The final results have been expressed in 
graphs which were used afterwards to give a 
graphical interpretation of the kinetic measure- 
ments, a method not previously used in this field. 

The object of the investigation to be discussed 
in this series of articles is to describe the behaviour 
of ion-exchangers in the preparation of boiler- 
feed water. As the treatment has been kept as 
general as possible the conclusions are also 
applicable to other fields. 


2. RELATIVE DATA FROM THE THEORY 


or CHROMATOGRAPHY 


To obtain a clear insight into the phenomena 
taking place in a column filled with ion-exchanger, 
it is useful to fall back upon those considerations 
in the theory of chromatography (de Vautrt [13], 
Guiueckaur [7]) in which these phenomena are 
brought into connection with the relative adsorp- 
tion isotherm(s). Starting from the well-known 
differential material balance (for piston-flow) : 


bc bc\ 
+7(*) =0 (1) 
Sx), bt), bt), 
it can be derived that in the case of instantaneous 
establishment of equilibrium in the column, the 


rate of displacement v, of a given concentration 
ec satisfies the relation : 


(=) : (2) 
bt Y + Sf" (ec) 

Here f‘(c) is the derivative of the isotherm 
q =f (c). Differentiation of equation 2 gives : 


bar 
de iy +f’ (c)}* 


This equation is very useful for describing the 
shape, or changes in the shape, of a chromatogram. 
When the isotherm is linear the rate of displace- 
ment v, is independent of concentration ; a given 
chromatogram retains its shape. 

When the curvature of the isotherm is positive 
(Fig. 1), loading of an empty column will result 
in a diffuse front since v, decreases with increasing 


Tangents and chords, 


Fig. 2. f’ (c) < 0. Tangents and chords, 

concentration. If, on the other hand, the isotherm 
shows a negative curvature (Fig. 2), v, would, 
according to equation 8, increase with increasing 
concentration which would lead to an impossible 
situation. It must be concluded, therefore, that 
the front is sharp and that all concentrations 


move at the same rate: 


(4) 


When the column is eluted with a pure solvent, 
an isotherm with a positive curvature will, in 
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Fie. 3 


— 


(b) 


Development of chromatogram for the case f” (c) > 0. 


(b) 


Development of chromatogram for the case f” (c) 


analogy with the foregoing, deve lop a sharp rear 


boundary, whereas an isotherm with a negative 


curvature will develop a diffuse one (Figs. 3a 
and 4a). 

With either type of isotherm the rear catches 
the the has a 
positive curvature the rate of c, in the (sharp) 
(diffuse) 
through 


up with front: when isotherm 


rear is higher than that of c, in the 
front the the chord 
(co. f (€9)) and (0,0) is smaller than that of the 


since slope of 
tangent in (¢y, f(¢,) to the isotherm see equations 
(2) (4) Fig. 1). Mutatis 
this also holds for the isotherms with negative 
The further concentration 
changes are shown in Fig. (3)b, (c) and 4(b), (c). 

If the isotherm shows an inflection point as in 
Fig. 5 > 90), 
equilibrium proceeds differently. 
concluded from the foregoing that with ¢ < ¢, 
the front tends to adopt a sharp shape and with 


and and mutandis, 


curvatures (Fig. 2). 


loading at instantaneous 


It might be 


Fic. 5. Isotherm with f (c) > 0 0 if ¢ ). 
c >, a diffuse one. However, the slope of the 
tangent (f’ (¢,)) in ¢, is smaller than that of the 
chord for (0-0) and f (¢,)); 


rate of ¢, in the sharp part of the front would be 


as a result the 


lower than that in the diffuse part. Consequently, 
the front will be straightened up to ¢,, which 
point is defined by : 
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(b) 


Fic. 6. 


Here again elution leads to rear formation as 
described above and finally the chromatogram is 
transformed as shown in Fig. 6(a), (b), (c). The rear 
the 
When the two sharp parts have caught up with 


overtakes front and ¢,, increases (c,,’). 


each other, c,, is given by the tangent drawn 
from ¢c,, to the top of the isotherm : 


When this point has been reached the further 
course of the curve is as shown in Figs. 6b and ce. 


Development of chromatogram for the case f’”’ (c) 


Fie, 7. Isotherm with (e) 0 if = ep). 


(b) 


Fic. 8. 


An analogous description can be given for the 
ease f'’’ (c) << 0 (Fig. 7). In Fig. 8 (a,b,c) the 
development of the chromatogram is indicated. 


Ion-EXCHANGE AS A FoRM OF 
ADSORPTIVE PERCOLATION 


some 
An 
example of ion-exchange reactions is the equili- 


lon-exchange in a column bears 


resemblance to a chromatographic process. 


Development of chromatogram for the case yr" tc) < @. 


brium between a sodium chloride solution and a 
cation-exchanger in the H*— form : 

NaCl + HR = NaR + HCl (6) 
or Na; Hy, == Na; + (7) 


In many theories these equilibria are represented 
by a formula ; however, these will not be discussed 
here because, like the empirical formulations, they 
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have only a limited validity. Therefore, we shall 
restrict ourselves to giving a graphical representa- 
tion of the equilibria in a following paper. With 
the above equilibrium the loading of the ion- 
exchanger q,), is a function of the concentrations 
of Na* and H* in the solution : 


=S ) (8) 


During the ion-exchange process the total 
concentration of the solution remains constant 
since every Na*-ion leaving the solution is 
replaced by an H °-ion. 

This is why 9, for different total concentra- 
tions in solution (a) is represented as a function 
of the equivalent fraction* of one ion type 
(e.g. Na*) in solution (<) : 


Ga = Sa (9) 


The formulae and observations mentioned above 
then continue to be valid provided ¢ is replaced by 
2a and it is borne in mine that in a column 
experiment a is constant. The material balance 
(equation 1) then changes into : 


av ay 8 (10) 

the equation for the rate of displacement of 
the equivalent fraction in a diffuse part of the 
chromatogram at instantaneous equilibrium 
(equation 2) into : 


ba ar (11) 
ay+Sf, (2) 

and the corresponding equation for a sharp 
part (equation 4) into: 


*The equivalent fraction in ion-mixtures of equal 
valency is of course equal to the molecular fraction 
this does not apply, however, when the ions differ in 


valency. 


Apam B. A. and Houmes E. L. J. So 
p’ G. U.S. Pat. 2,366,008 


Bapvour R. F. and E 
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dx av 


dt 4 


(12) 


Here a, = the equivalent fraction of Na” -ions 
in the liquid which is flowing in. 

In analogy to the substitution of aa for e¢, 
q, might be replaced by the equivalent fraction 
of the ion in the exchanger multiplied by the total 
capacity ; as the latter quantity is often difficult 
to determine this substitution will not be made 
here. 

With these equations, the ion-exchange 
mechanism has been expressed in formulae 
which are normally used in chromatography. 


Acknowledgement—The authors thank Dr. C. van 
Hererpven for valuable discussions about this 


subject. 
NOTATION 
a total ion concentration in solution (M L-4) 
‘ concentration in the liquid (M L-) 
fie) concentration in the absorbent at 
equilibrium (M L-) 
Sq (2) concentration in the ion-exchanger 
at equilibrium (ML) 
q = concentration in the absorbent (M L-) 
Yq = concentration in the ion-exchanger (M L-) 
Ss cross-sectional area of the column (L?) 
t time (T) 
t rate of flow (L3 T=) 
v, = displacement rate of the concentration ¢ 
in the column (LT) 
r distance to the point where the liquid 
enters the column (L) 


a = equivalent fraction of one ion type in 
solution 
free space per unit length in the column (L?) 


Subscripts 
0 condition before the liquid enters the column 
b inflection point in the isotherm 
L = liquid 
R = ion-exchanger 
r tangent point 
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Abstract—A description is given of current methods for the determination of the equilibria 
of ion exchangers in electrolyte solutions, by means of which the behaviour of a column can be 
accounted for ; among these, the “ column method” is to be preferred. However, it is found 
that experiments carried out according to the column method easily lead to results that are 
inaccurate or even wrong. It is shown that this disadvantage can be eliminated by the use 
of buffer solutions having the ion concentration for which the equilibrium is required, for a buffering 
anion has a great influence on the way in which the equilibrium is reached with cation exchangers, 


and the same influence is exerted by a buffering cation on anion exchangers. 


The determination of a number of equilibria with cation and anion exchangers is described 
with reference to the purification of water. It is striking that in the resulting isotherms inflection 
points occur very frequently ; these points are attributable to the polyfun tionality of the 
products. In general, an isotherm may have a great variety of shapes, dependent on the structure 
of the exchanger and on the composition of the solution. A classification of isotherms on this 
basis is given. The exchange of uni — univalent ions proves to be independent of the total ion 
concentration of the solution. 


Résumé—On donne une description de méthodes courantes pour la détermination des équilibres 
d’échangeurs d’ions dans des solutions électrolytiques, équilibres qui permettent de calculer 
le comportement d'une colonne ; de toutes ces méthodes on préfére appliquer la méthode de 


la colonne. Cependant, on a constaté que les expérimentations effectuées selon la méthode de 


la colonne entrainent facilement des résultats inexacts et méme faux. 


On démontre que ce désavantage peut-étre évité par l'emploi de solutions tampons qui ont 
une concentration d’ions pour laquelle on désire connaitre l'équilibre. C'est qu'un anion tampon 


exerce une grande influence sur la maniére dont on atteint léquilibre avee des échangeurs de 


cations, et un cation tampon exerce la méme influence sur les échangeurs d’anions. 


La détermination d'un nombre d’équilibres avec des échangeurs de cations et d’anions se 


trouve décrite par rapport a la purification d'eau. Il est remarquable que les isothermes qui en 


résultent manifestent trés souvent des points d’inflexion, points qui sont dus 4 la polyfonctionnalité 


des produits. En général, une isotherme peut avoir une grande diversité de formes en fonction 


de la structure de léchangeur et de la composition de la solution. On donne une classification 


disothermes a cette base. 


L’échange d’ions uni - univalents s'avére indépendante de la concentration totale d’ions de la 


solution. 


Zusammenfassung—bDie heute iblichen Methoden fiir die Bestimmung der Gleichgewichte von 
Jonenaustauschern in elektrolytischen Losungen werden beschrieben, mit deren Hilfe das Verhalten 


einer Kolonne ermittelt werden kann. Unter diesen wird die * Sdéulenmethode ” bevorzugt. Man 


findet indessen, dass Versuche, die gemiiss der Siulenmethode ausgefiihrt werden, leicht zu unge- 


nauen oder sogar falschen LErgebnissen fiihren. Es zeigt sich, dass dieser Nachteil durch die 


Verwendung von Pufferldsungen vermieden werden kann, welche diejenige lonenkonzentration 
J 


besitzen, fiir die das Gleichgewicht gesucht wird. Ein pufferndes Anion hat einen grossen Einfluss 


auf die Art, in der das Gleichgewicht mit Kationenaustauschern erreicht wird, und der gleiche 


Einfluss wird durch ein pufferndes Kation auf Anionaustauscher ausgeiibt. 
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Die Verfasser beschreiben die Bestimmung einer Anzahl von Gleichgewichten mit Kationen- 


und Anionen-Austauschern im Hinblick auf die Wasserreinigung. 


Es fallt auf, dass die sich 


ergebenden Isothermen sehr hiufig Wendepunkte besitzen. Diese Punkte lassen sich durch die 


sehr verschiedene Natur der Produkte erkliren. 


Im allgemeinen kann eine Isotherme sehr 


verschiedene Formen haben, abhiingig von der Struktur des Austauschers und der Zusammen- 


setzung der Lésung. Auf dieser Grundlage wird eine Klassifikation von Isothermen mitgeteilt. 


Der Austausch von ein-einwertigen Ionen erweist sich als unabhingig von der gesamten Ionen- 


konzentration der Lésung. 


In part I it has been proved that in chromato- 
graphy, in the case of an instantaneous equilibrium 
in the column, the change in the shape of the 
chromatogram can be quantitatively and qualita- 
tively predicted from the shape of the adsorption 
isotherm. A similar observation applies to ion 
exchange. 

In reality, the establishment of the equili- 
brium is, of course, not instantaneous, but even 
so, knowledge of the position of the equilibrium 
is of the utmost importance for calculating the 
driving force of the process. 

In this part, attention will be paid to the 
determination of the equilibria and to the classi- 
fication of ion exchangers according to the shapes 
of the equilibrium curves. 


METHODS OF DETERMINING 
EQUILIBRIUM 
The relation ¢, = f,(«) can be determined by 
the following methods : 
(1.1) The batch method 
This method was used for instance by Grecor 


and Breoman [8]. While 
tinuously, a quantity of 


being stirred con- 


given exchanger is 
allowed to reach equilibrium with a_ solution 
containing the ions of which the equilibrium is 
to be determined. Once equilibrium has been 
reached, the compositions of the solution and 
the exchanger are determined. The disadvantage 
of this method is that it is impossible to predict 
the value of « which the solution will eventually 
have, so that it is often necessary to carry out 
several experiments in order to scan a certain 
range. 


(1.2) The derivation from a continuous elution 
curve 

This 
elaborated by Giurckaur [5], [6], Sinuén [15], 


less widely used method has been 


[16] and Yanowsk1 [18]. By integration of the 
following equation (equation 11 in 
paper [10}) : 


our first 


ai 


1) 

ay + Sf, (a) ( 
a relation is obtained between f.(«) and the 
volume of liquid passed through the column. 
From the elution curve « can be determined as 
a function of this volume, and by graphical 


integration of f, (x) it is then possible to deter- 


mine f, (x). 
The this 
pointed out by GLueckaur and Coares [7], is 


disadvantage of method, already 
the probability of rather serious errors in the case 
little 


If for some reason the chromatogram is disturbed, 


of isotherms deviating from linearity. 
or if a kinetic process is involved, such a distur- 
bance may give rise to the erroneous conclusion 
that there is a diffuse front, as with linear iso- 


therms all concentrations have the same velocity. 


(1.38) The column method 

A solution having a known value of « is passed 
over a column containing an ion exchanger until 
the composition of the effluent becomes equal 
to that of the incoming solution. A disadvantage 
is that, if a widening front is involved (with a 
sharp front the experiment in general proceeds 
rapidly), there is a chance that a large amount of 
liquid has to be passed over before the correct 
equilibrium is reached; for, if a certain « is 
accompanied by a high y a (x), the rate at which x 
the 
to equation (1), be low. 


moves through column will, according 


That this is not an imaginary danger may be 
seen from the fact that, whereas the Na” H* 
equilibrium for a monofunctional strongly acid 
exchanger can easily be determined with a 
solution of hydrochloric acid and common salt, 
this method cannot be used for a polyfunctional 
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exchanger possessing both strongly acid and 
weakly acid groups. Later it will be seen that 
with the latter type of exchangers i (x) is high 
if « has a high value. One cause is that, on the 
strength of analysis, there is a tendency to think 
that the effluent has the same composition as the 
incoming solution, while owing to a small error 
in the analysis a large error in f, (x) may occur if 
fa (a) has a high value. 

It will be that the latter 
exchangers the equilibria may well be determined, 


shown also with 


provided use is made of buffer solutions. 


(1.4) The column method using buffer solutions 


In itself the use of buffer solutions is not a 
novel technique : these solutions have long been 
used in the determination of the relation between 
pH and total capacity [8]. In the 
buffer are to be 


determining equilibrium curves, so that care 


present case, 


however, solutions used for 
will have to be taken that the « values show a 
wide spread and that the total ion concentration 
in the inflowing solution in the various tests is 


always the same and equal to the ion concentra- 


tion for which the equilibrium curve is to be 
determimed. 
However, if buffer solutions are used, not 
only the equilibrium : 
Na, Na, (2) 


but also the equilibrium : 


plays a part, which implies that in the column 


Cre Cn is no longer’ constant, but 
Na CH 
In addition to the previously determined equiva- 


lent fraction : 
(4) 


in which the index 1 relates to the isotherm to 
the auxiliary quantity 


° 
= (C wat) wat TC pt) 


be determined (Fig. 1), 
Xe is now introduced : 
(5) 


corresponding to a Se (%) 
(ef. Fig. 1, curve 2); this indicates by its shape 


(Cxa+)/ (nat + Cyt Cyx) 


second isotherm 
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Fic. 1. Curves of types 1 and 2. 


whether or not a sharp front is formed in the 
column. 

A sharp front (/..' (a) < 0) is ideal; the only 
requirement is that f..(%,) be not too great ; in 
any case x, must be small, and hence, in accor- 
dance with equation (5), ¢y, great. This re- 
quirement can be fulfilled by a correct choice 
of HX 
ently weak. 

It will be clear that without a buffer f, («,°) 
(x,°) ts very great, 


i.e. by ensuring that this acid is suffici- 


cannot be reached, because f,' 
whereas, if a buffer is used, this load can be 
realized (f,(a,°) = fy (a_°)), because fy ts 
small. Each point to be determined of the iso- 
therm 1, i.e. each value of the equivalent fraction 
in the solution («,°), has its own curve 2, whose 
shape indicates how the f,(x) value correspond- 
ing to «,° is reached; consequently there is a 
series of curves of type 2 relating to each iso- 


therm to be determined. 


The above observations can be illustrated by describing 
how a second isotherm can be calculated if the first should 


be known. For equilibrium (3) the following equation 


holds : 
K = ¢x~)/(CHx) (6) 
and, owing to the electric neutrality of the solution : 
+ Cyt = (7) 
Finally it is assumed that : 
+ Cut + CHx = (8) 


If ex ,+ and cyyy are eliminated from equation (5) by means 
of equations (6), (7) and (8), it follows that : 


| 
vO 
1957 
| 


K — “nt (9) 


Next, cyy+ is chosen and a, calculated from equation (9), 
bearing in mind that to ag there belongs the same f(z) 
as to: 

ay = c°)/(ag c° + (10) 


It is clear that ag is smaller according as the dissociation 
constant K is lower, and depending on the value of K, 
calculation results in an isotherm similar in shape to curve 
2 in Fig. 1. 

An additional advantage of buffer solutions is 
that «,° is more strongly fixed, which contributes 
to the accuracy of the experiment, particularly 
for high « (low H* ion concentration) and f’ (x) ; 
this advantage is also obtained if the batch 
method is used. 

From these considerations if is concluded that 
the column method with buffer solutions is to 
be preferred to all others. 


2. EXPERIMENTS TO DETERMINE 
THE EQUILIBRIUM 

(2.1) Cation exchangers 

Of the cations 

Ca** and Na* are the most numerous and qualita- 


contained in natural waters 


Studies on ion-exchange—II 


No’ conc.0°003 N 
Ca”*conc. 0-01 N 


Ca’*;conc. 0-003 N | 


the highest and Na~ the lowest affinity for cation 
exchangers of all the positive ions contained in 
water [12]. Equilibria determined for 
these two kinds of ions. The total ion concentra- 


were 


tions chosen were 0-01 and 0-003 g equiv/l; these 

concentrations represent about the maximum and 

mean salt contents of natural waters. The follow- 

ing exchangers were examined [2]. 

(a) Imac C 12, 
corresponding to Dowex 50 and Amberlite 
IR 120. 
strongly acid exchanger. 

(b) Dusarit 
strongly acid and weakly acid groups. 


a sulphonated polystyrene, 


It is considered a monofunctional 


a sulphonated coal with both 


(c) Dusarit VKG, likewise a sulphonated coal, 
which has, however, more weakly acid and 
fewer strongly acid groups. 

(d) Lewatit, a polyfunctional exchanger on a 
synthetic resin basis. 

(e) Stamion, a weakly acid exchanger obtained 
by a coal with a 


treatment of oxidized 


solution of sodium hydroxide [4]. 

(f) Amberlite IRC 50, a monofunctional weakly 
acid exchanger, probably a copolymer of 
metacrylic acid and divinyl benzene [3]. 


Samples consisting 


of 25-35 ml of ion ex- 
changer were treated 
in small columns con- 
solutions of 
buffer salts, until the 
effluent had the same 
the 


taining 


composition as 


incoming solution, 
| which was checked 


by measuring the pH 
value. 
The buffer solutions 
were prepared as 
If the equili- 
brium curve for Na* 


follows. 


ions was to be deter- 
mined at a total con- 


tively the most important, as usually Ca** has 
2000, 2000 
~ 
E 
1500) 
| 
4 | 
3 3 
No’ iconc.0-O03N 
Co”conc.0-O1N 
o—-—o Ca**iconc. 0-003 N 
re) 
02 04 O06 O8 1-0 2 


Fic. 2a. Equilibrium curves of Imac C 12. 


Fic. 2b. Capacity-pH curves of Imac C 12. 
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centration of Na* and 
ions = 0-003 N, 
(where N_ (normal) 


4 § 6 7 8 
pH 
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7 
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* 
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e——e Nod ;cone 


-) 
Fic. 38a. Equilibrium curves of Dusarit 5. 


refers to lg equiy l.). use was made, (e.g. for 
a value of « 0-6), of a solution of sodium 
0-0018 N. By add- 
ing 1 : 1 phosphoric acid the solution was brought 
to a pH log 0.0038 (1 0-6) 2-92. The 


adjustment was made by means of a plHl-meter 


hydroxide of 0-6 « 0-008 


(radiometer). For the concentration 0-01 N an 
analogous method of preparation may be used. 
The Ca®* buffers were prepared from caletum 
acetate solutions. 

After the experiment the exchanger was 


regenerated with excess hydrochloric acid; the 


Fic. 4a. Equilibrium curves of Dusarit VKG. 
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pH 


Fic. 3b. Capacity-pH curves of Dusarit 5S. 


Na’ ion in the regenerating solutions was deter- 
mined by precipitation with zincuranylacctate, 
and the Ca* ion by titration, for which the 
complexone method was used. The f, value 
corresponding to the « value given above could 
then be calculated directly as m equiv 

The results are given in Figs. 2a—7a which show 
the loads of the exchangers as functions of the 
equivalent fractions ay, and a, ; as these graphs 
are not very clear at high values of «, in Figs. 


2b-7b the loads are shown as functions of the pH 
value. 


Fic. 4b. Capacity-pH curves of Dusarit VKG. 


4 ©-003N z e——e No”;conc. 0-003 N 
E E Ca’*:conc. O-OIN 
1000} 1000) 
y 
j 
If 
\ 0 
nee E J 
= 0 ty 1000} 
= A : 
/ a / 0-003 N 
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am Ca** conc. O-C 


Fic. 5a. Equilibrium curves of Lewatit. 


06 O08 


a 


Fic. 6a. Equilibrium curves of Stamion. 


(2.2) Anion exchangers 
The 


exchangers is, in principle, subject to the same 


measurement of equilibria for anion 


considerations as those given for cation 


exchangers; in this case also buffer solutions 


have to be used. There are, however, two 


practical difficulties. 
In the first place it is the cation that has to 


f 
2000) 


| 
1500}- 


Fic. 5b, Capacity-pH curves of Lewatit. 


1500,—__— 


No’ :conc 


2 


Fic. 6b. Capacity-pH curves of Stamion. 


act as a buffer with anion exchangers, and the 
choice of cations is limited. After some searching 
we found as suitable cations piperidine and 
homo-piperidine, the basicity of which corresponds 
to the acidity of phosphoric acid, and mono- 
ethanolamine with a basicity corresponding to 
the acidity of acetic acid. 


In the second place the pH adjustment of the 
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solutions to be used for determining the equilibria 

in the alkaline range cannot be made by means of 

a glass electrode, owing to the wide variation 

in the range covered. With a hydrogen electrode 

reproducible results were obtained. The OH Cl 
equilibria were measured, and those for the 
following anion exchangers are mentioned [2] : 

(a) Amberlite IRA 401, an exchanger with 
quaternary (alkylated) ammonium groups 
on a polystyrene basis. 

(b) Permutit ES, like (a), but not monofunctional 
as, in addition to alkyl groups, there are 
probably also alkylol groups (corresponding 
to Dowex 2 and Amberlite IRA 410), 


(ec) Dowex 3, like (a), but containing only 
weakly basic groups. 

The equilibria were determined in the same way 
as for cation exchangers. Solutions with total 
OH™ and Cl” ion concentrations of 0-01 N, were 
prepared by adjusting solutions of « 0-O1N 
hydrochloric acid to a pH value of 14-log 0-01 
(1 a) by addition of monoethanolamine (for 
lower OH” ion concentrations) and piperidine or 
homo-piperidine (for higher OH ion con- 
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Equilibrium curves of Amberlite IRC 50. 
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o—-—oCa’ “conc. 0-OO3N 
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Fic. 7b. 
Capacity-pH curves of Amberlite IRC 50. 


centrations). After regeneration with an excess 
amount of sodium hydroxide solution the Cl~ ion 
concentration was determined by titration accord- 
ing to Mohr. Results are shown in Figs. 8 and 9. 

Not so much for the characterization as for 
the practice of ion exchange it is desirable to 
know also the acid adsorption equilibria for 
anion exchangers. Some experiments which 
were made to determine the hydrochloric acid 
equilibria for Dowex 3 and Permutit ES will 
therefore be described. 

The determination of these equilibria was car- 
ried out in a similar manner to that used for the 
equilibria, except that instead of 
buffer solutions, solutions of hydrochloric acid 
(which may be considered as limiting cases of 
buffer solutions) of various concentrations were 
used. 

It was found that with Dowex 3 the equilibrium 
was only established very slowly: even after 
some days had elapsed, the concentration of the 
hydrochloric acid in the effluent was lower than 
that in the incoming solution. Therefore the 
method was slightly altered ; first the exchanger 
was loaded with 1N hydrochloric acid and 


K. Kramer, J. C. 
' ~ 
4 No’ cone. O-OTN 
E 
e——e ;conc. 0-O03N 
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Fic. 8. 
Equilibrium curves of anion exchangers. 


next with solutions of lower concentrations 
(< 0-02 N) for which the isotherm had to be 
determined. With the altered method the con- 
centration in the effluent (which was at first 
higher) after a reasonable time became equal 


to the concentration in the fresh solution. An 


Dowex 3 


o—ePermutit &S 


Fic. 10. 
Acid adsorption curves of Dowex 3 and Permutit ES. 


———.» Permutit ES 


e—-—eDowex 3 


Fic. 9. 
Capacity-pH curves of anion exchangers. 


explanation of the effect of this alteration of 
the method of determining the equilibrium will 
be given at a later stage. 

The hydrochloric acid adsorption equilibria for 
Dowex 3 and Permutit ES are given in Fig. 10, 
(2.3) Discussion 


(a) 


shows a 


Inspection of the isotherms of Figs. 2-8 
of inflection 
GLUECKAUF was the first to find this type of 
for exchangers, viz. with the 
equilibrium Cu*?* — H®* for Zeo-Karb which is 
presumably a polyfunctional exchanger [6]. Fur 
ther examples are the equilibria Cl NO, for 
Dowex 2, found by Wuearon and Bauman [17] 
(cf. Zeecers [19)). that 
this type of curve will be found where the affinities 


striking number points. 


isotherm ion 


GLUECKAUF suggests 
of the competing ions for the exchanger are of the 
same order. By way of example he refers to the 
systems examined by Rorumunp and 
[13], viz. Ag* — Tl* and Na NH, for alumin 
ium silicate exchangers. Against GLUECKAUF’s 


opinion it may be adduced that, if the affinities 


of the various ions for the exchanger are of the 
same order, a linear isotherm is to be expected. 
If the experiments of Roruamunp and Kornre_p 
are plotted graphically it is, indeed, seen that 
linear isotherms in the « 


these authors found 


region they investigated. 


211 


1500 
E 
WA 
1000 — 1000} A + | 
= | | 
| 
a DH 
VOL. 
057 / 52 
~ | 
| | 
° 4 6 2 6 20 
c megq/t 
= 


K. Kuamer, J. C. H. 


Like SELKE [14] the authors believe that this type 
of curve is inherently bound up u ith the poly- 
functionality of an exchanger. 


From Figs. 2-7 the impression is obtained that 
the curves for Dusarit VKG, for example, might 
have been formed by superposition of the curves 
for Imac C 12 on those for Stamion. A mixture of 
exchangers such as Imac C 12 and Stamion, or 
an exchanger containing acid groups of different 
strengths will, in the case of equilibria in which 
H* ions take part, give curves with inflection 
points. It may be said that if the metal ion in 
solution has a small equivalent fraction, when 
the strongly acid groups on the exchanger ar 
chiefly involved, the affinity of the metal ion 
for the exchanger is smaller than that of the H- 
ion, whereas at high values of «, when the weakly 
acid groups also play a part, the aflinity of the 
metal ion for the exchanger will become pre domin- 
ant. The course of the equilibrium curves of Imac 
C 12 at high «, and the analogy with the poly- 
that 


number of 


functional exchangers, suggest this ex- 


contains a_ certain 
This 


agree with what Boyp [1] presumed in connec- 


changer also 


weakly acid groups. conclusion would 
tion with the sulphonated polystyrenes. 

(b) Mathematical representation of equilibrium 
In general it is impossible to represent 


The only 


curves. 
the equilibrium curves in a formula. 
regularity is that the Na* H 


in most cases independent of the total concen- 


equilibrium is 


tration in the solution and hence exclusively 
depends on the proportion between the ion 
concentrations. 

This has been observed by Honpa [9] among 
others. Perhaps this relation has a deeper cause. 
If we 
equilibrium, the following equilibria must be 


assume the occurrence of a Donnan 


established in a system containing Na’, i 


and Cl 


1OnS 


* ¢ Cut * Caz 
Hp Hp “ay (11) 
Cap * cl (12) 


L 


from which follows, after division : 


(13) 


Cup CNa} CHy/ 


If it is further assumed that cup depends only 
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on the load of Na* ions (Cra), it will immediately 
be clear that : 


naj, = S (exat/euz) (14) 


Unfortunately, Cup can only be calculated from 
equation (13), and cannot be determined indepen- 
dently, as the number of non-dissociated groups 
on the exchanger is not known. 

If equation (13) is assumed to be correct one 
may get an impression of the dissociation constants 
of the acid groups on the exchanger in dependence 
On from and from a determination of 
the total number of exchangeable groups q, 


using the equation : 


K = (CNaj, 


(naj, — (15) 


| 


Table 1. 
Dussarit VAG 


K K 
(g mol/L.) (g mol/L.) 
630 
0-99 760 
900 
1070 


0-092 
O-O13 
0-0027 
0-00072 


0-999 


For 


corresponding to equation (14), holds : 


anion exchangers the following relation, 


i, Comp) 


and, for a preponderantly acid medium : 


f (16) 


which corresponds to the results of Kunin and 
Myers [11] and to the results of the present 
authors (Fig. 10). 

(c) The position of the capacity-pH curves 


depends on the total ion concentration in solution. 


Grecor and BrecMan [8] came to the above 
conclusion, and therefore used high salt concen- 
trations in solution, with the result that the capa- 
city-pH curves became relatively less dependent 
onthe salt content. This effect is easily understood; 
from Fig. 2b-7b it can be seen that at a high pH 
value the capacity usually becomes less dependent 
on the salt concentration. If now the concen- 
tration is raised, the capacity-pH curve shifts 
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— Types of ion exchanger: 

Reacting solutions 


cotion exchanger 


nd 


Neutral salts 


Buffer solu tions 
janion buffer [ABs] 
lcation buffer 


Acids and bases 


Reaction of fresh ion exchanger with 


(caustic solution 


locia solution [aS) 


Regeneration of 
oaded 


mpietely 


(acid solut 


coustic solutior 


Fic. 11 
to the left. Doubling the concentration results 
in a shift of the curve of 0-3 pH unit ( 
so that a measurement in the normal pH range 
(pH > 2) then 


relatively little dependent on the total concen- 


log 2), 


will show the capacity to be 
tration in solution. 

The use of high salt concentrations provides 
comparative information on the acidity distribution 
of the groups in the exchanger. This information 
cannot be used for calculating equilibrium curves, 
as it gives capacity values for very high « only. 

3. CLASSIFI<C Ion 


AND OF 


ATION OF 


EXCHANGERS Ion EXCHANGE 


ISOTHERMS 


On the basis of the foregoing it is possible to 
classify ion exchangers according to the shapes 
of their equilibrium curves ; this classification is 
shown in Fig. 11. 

It is seen that, depending on the nature of the 
exchanger (monofunctional, strong or weak ; or 
polyfunctional) and of the type of the solution 
(neutral salts, acids and bases, buffered or 
unbuffered), and on the question whether loading 
or regeneration takes place, a great variety of 


isotherms can be obtained. 


lon exchanger 


Polyfunctional Monofunctional 
ication exchanger 


janron exchanger 


{acid cation exchanger 
)basic anion exchanger 


Weakly 


and 


ABS) 
ond 
prac 


1 


and 
prac las] 


Classification of ion exchange equilibrium curves. 


In Part III it will be seen what quantitative 
data on the capacity follow immediately from 
the shape of the equilibrium curve. 


NOTATION 


(ML-) 


total ion concentration in solution 
concentration in the liquid 
( for definition, see eq. (8) 


f(x) 


concentration in ion exchanger at 
(ML-%) 
(ML-) 


equilibrium 
kK dissociation constant 


total number of groups in the exchanger 
(ML) 


(L?) 
t time (T) 
T-) 


(whether or not dissociated) 


cross-sectional area of column 


v = rate of 


x = distance to the point where the liquid 
enters the column (L) 


equivalent fraction of a given ion 
a 


in solution 


free space per unit length in the column (L?) 
Indices and subscripts 

L = liquid 

R = exchanger 


7) condition before the liquid enters the column 


—— 
: 
(x) (a | f (ce) 
| 4 saCE|nS 4 jpfCE ins) | 
| | | and ond | end | 
| + wy or 
| SbAE|nS prAE|ns | | 
— 
| and |! and 
| SDAE|CBS | | 
= 
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The hydrodynamics of flow between horizontal concentric cylinders—I 
Flow due to rotation of cylinder 
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Abstract 
centric cylinders with the inner cylinder only rotating have been studied for liquid levels other 


The approximate velocity distributions for laminar flow between horizontal con- 


than when the system is completely full. The stability of the laminar flow has been shown to 
depend on the parameters S and dy/d as used by Di Prima, where dp is the distance of the 
surface of zero velocity from the inner cylinder and d is the width of the gap. The secondary 
motion resulting from instability has been investigated and the wavelength of the vortices 
produced between the inner cylinder and dp agrees with that for the same value of dy/d due to 
rotation of both cylinders in opposite directions when the annulus is full of liquid. 


Résumé 
écoulement laminaire entre des cylindres concentriques horizontaux, le cylindre intérieur 


Les auteurs étudient la répartition approximative des vitesses dans le cas d'un 


tournant. Les mesures ont été faites pour des niveaux différents de liquide, le systéme n’étant 
jamais enti¢rement noyé. Ils ont montré que la stabilité de cet écoulement laminaire dépend des 
paramétres S et d)/d, comme l'a indiqué Di Prima: d) est la distance de la surface de vitesse 
nulle, le cylindre intérieur étant pris comme origine, et d la largeur de espace annulaire. 

Le mouvement secondaire résultant de linstabilité a été étudié, et la longueur d’onde du 
tourbillon produit entre le cylindre intérieur et dy, concorde avec celle correspondant & la méme 
valeur de d,/d, dans le cas de la rotation de deux cylindres en sens inverse, quand l’éspace 
annulaire est remplie de liquide. 


Zusammenfassung—Es wurden die angeniherten Geschwindigkeitsverteilungen fiir laminare 
Stromung zwischen waagerechten, konzentrischen Zylindern, von denen nur der innere rotierte, 
bei nicht vollstindig gefiilltem System untersucht. Die Stabilitat der laminaren Stromung erwies 
sich als abhingig von den Parametern S und d,/d, wie von Di Prima angegeben, wobei dy 
die Entfernung der Schicht mit der Geschwindigkeit Null vom inneren Zylinder und d die 
Spaltweite ist. Die von der Instabilitét herriihrende Sekundiirbewegung wurde untersucht ; 
die Wellenliinge der zwischen dem Innenzylinder und d, erzeugten Wirbel stimmte mit jener fiir 
gleiche d,/d-Werte bei entgegengesetzter Rotation beider Zylinder und vollstindig gefiilltem 
Ringspalt tberein. 


INTRODUCTION making, a sheet of fibre is taken off a suspension 


Tue flow of liquids in curved passages is of on to a drum covered with wire gauze which 
rotates in a vat full of a fibre suspension. 


interest in a number of chemical engineering 
problems. A rotating drum filter, for instance, 
consists of a large cylinder which rotates inside 
a container, the suspension to be filtered being 
in the annulus. 


plicated form of motion : it is rotated by the drum, 


The suspension has a com- 


is pumped through the annulus under a positive 
head, and is finally given a radial component in 
the filtering region. Similarly, in paper- and board- 


In these and similar flow systems there are 
(1) the distribution of 
velocity and pressure heads across and along the 


essentially two problems : 


system under different regimes of flow; and (2) 
the critical speeds at which the flow changes from 
laminar to a more complicated pattern. For 
example, in paper-making, longitudinal streaks 
‘an arise in the formed sheet when a certain type 
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of vortex motion predominates in the flow, due 


to its instability in a curved passage. 

In this paper will be discussed the motion of a 
liquid contained between two concentric cylinders 
with the inner cylinder rotating. The axis of the 
cylinders is horizontal and a free liquid surface is 
present. 


THEORY oF FLOW BETWEEN ROTATING 


CONCENTRIC 


CYLINDERS 


It will be assumed that the annulus between two 
concentric cylinders contains an incompressible, 
homogeneous Newtonian fluid. The problem 
consists of determining the character of the 
laminar motion existing at slow speeds and 
secondly the speed at which this motion becomes 
unstable. Consideration of the type of motion 
produced when instability occurs—vortex or 
turbulent motion—-is important but will only be 
briefly dealt with in this paper. 


Laminar Motion (slow speeds) 


When the annulus is completely full of liquid 
the mode of flow is well known and the tangential 
(Couette) velocity distribution between the 
cylinders is given by Lamp [1]. If the cylinders 
rotate in opposite directions there is a cylindrical 
surface of radius r, between the two cylinders at 


which the tangential velocity is zero. This is 


given by 


(1) 


If the common axis of the cylinders is horizontal 
it is possible to drain liquid from the annulus 
between the cylinders until the top of the inner 


“dry ”. The Couette velocity distribu- 


cylinder is 
tion will not longer hold between a and b for 
this system. 

Consider the flow when the outer cylinder is 
stationary and the inner rotates (yu 0). 

It will be assumed that as the rotating cylinder 
passes through the free horizontal liquid surface 
it carries with it a liquid film of negligible thick- 
ness compared with the annulus width d. The 
rotating liquid, next to the moving inner cylinder, 
on reaching the free horizontal liquid surface, 
cannot continue in the same direction and must 
return in the opposite direction as shown in Fig. 1. 
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Fig. 1. Rotation of inner cylinder with complete reversed 
flow. 


The free liquid surfaces at the two sides of the 
cylinder will be at different levels and Ah, the 


difference in height, will be a measure of the 
* reverse 


pressure head required to produce this 
flow.” 

This simple picture of a “ dry” top will only 
hold for liquids of low viscosity. If the viscosity 
is high a film of liquid of significant thickness will 
be carried over the top of the cylinder. It is 
hoped to deal with this more complicated case 
in later work. 

The velocity distribution between a and b will 


be approximately as shown by Fig. 2. 


| Ar 
- 


Fic. 2. Approximate velocity distribution between 
concentric cylinders with complete reversed flow. 


For values of d/a < a we may consider the 
problem as one of flow between parallel plates, 
one of which moves at velocity 2, a relative to 
the other. 

The velocity distribution [1] is given by 


2, a (2) 
( 2vp ade 
If the net flow is zero, 
d 
= }dQ,a — 4p 
12 vp ad@ 


\Y 
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2 vp ad@ d? 


When v = 0, from equations (2) and (3) : 


d,/d = 1/8 (4) 


These relationships apply when there is a 
complete reversal of flow. For a fluid of low 
viscosity it is possible to drain liquid from 
between cylinders until the liquid level lies 
between the highest points of the two cylinders 


and to obtain a steady state of flow. Rotation of 


the inner cylinder in this case will produce 
incomplete reverse flow as shown by Fig. 3*. 


Fic. 3, Rotation of inner cylinder with partially reversed 
flow. 


Assuming an approximate velocity distribution 
as in Fig. 4 and assuming that the velocity at 
d’ (Fig. 3) is the same for any point at this 
distance from the inner cylinder, the equation of 


| ! | 
| 
| 


b+ 


| 
| 
| 
| 
distribution between 


Fic. 4. Approximate velocity 
concentric cylinders with partially reversed flow. 


continuity is given for the region between d’ 
and d, in the body of the liquid and far from the 
free surface, by : 


*When the viscosity is high the * bubble” is pumped 
bodily or partly around between the cylinders. 
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Flow due to rotation of cylinder 


d, d 


| v dx = ] v dx 


d’ d, 


We shall assume from equation (4) that in this 


case 
l 2d’ 
or “o (1s) 4 (5) 
d 


Stability of Laminar Motion 
Full annulus 


The stability of laminar flow between rotating 
concentric cylinders has been studied by many 
workers. G. I. Tay.Lor [2] was the first to show 
analytically and experimentally that for cylinders 
of given radii and a liquid of a certain viscosity 
there were definite angular velocities of the outer 
and inner cylinders at which instability would 
occur. He further showed that the motion 
resulting from instability was three-dimensional 
and consisted of pairs of toroid-shaped vortices, 
the rotation within each being in the opposite 
direction to its neighbour. The spacing of these 
toroids along the axis of the cylinders is constant 
for a particular apparatus and speeds of rotation 
of the cylinders. The wavelength or width of 
adjacent pairs of vortices depends on the cylinder 
radii and speeds of rotation. 

TayLtor showed that when only the inner 
cylinder was rotating the critical angular velocity 
was given by 


Q2 a? 


0-0571 (6) 


0-00056 


where f = 1 — 0-652 (d/a) 


and d = (b — a). 


Lewis [3] confirmed TayLor’s work experi- 
mentally and showed that approximations involv- 
ing (d/a) < a used by Taytor in his theoretical 
analysis held for (d/a) as high as 0-7. 

The solution of the stability problem mathe- 
matically has been further studied by several 
workers. 

SyNnGE [4] has shown that for stability of a 
rotating fluid the circulation 


must increase 


| 
| 
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outwards. That is for angular velocities ,. 
w, at r, and r, the flow is stable if 


T 


When the cylinders rotate in opposite directions 
there is an inner section where the angular 
velocity is decreasing outwards and this section 
will tend towards instability. The outer section 
where the angular velocity is increasing outwards 
will tend towards stability. 


More recent work, mainly on the theory of 


instability is reviewed by Lin [5] and Orpennem 
and Hucues [6]. 

Dean [7] has shown theoretically that distur- 
bances of the same type as those described by 
TayLor [2] for Couette flow can exist for a 
parabolic velocity distribution. He showed that 
the critical velocity for instability is 


36 


There has been no direct experimental verification 
of this criterion, but it appears to explain why the 
Reynolds number for instability in a curved pipe 
is lower than that for a straight pipe. 

Di Prima [8] has plotted the results of expert- 
ments carried out by Tayior [2}. He has 
expressed the angular velocities (2, and 2,) in 
terms of 
822 [1 — (b*/a*)P de (7) 

and plotted experimental values of S against 
d,/d (see Fig. 9) where Q. is the critical angular 
velocity of the inner cylinder. 

We may express S in terms of d, by use of 
equation (1) 


8 122 


[4 + 8 (d, a) 5 (d,, (d, 


a)*| (7a) 


for dyja < 1 this may be approximated by 


s 2 QF a 


or Q (8) 
2a d, d, 


The wavelength of the vortex motion parallel 
to the axis of the cylinders resulting from 
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instability is given by Fig. 5 where ~ is the 


wavelength [8]. 


+ 
>} 
exer tol = 3-80 
x imern 
b= 4-035 
= — 


O86 08 TO 


qd, a 


Fic. 5. Wavelength of the vortex motion. 


For complete reverse flow, SyNGr’'s criterion 
[4) would suggest that the inner (Couette) 
section tends towards instability. The outer 
section, not having a simple Couette velocity 
distribution, but a parabolic distribution, may 
affect the stability of the inner section in a 
different way from Couette flow and may tend 
towards instability itself. 

In order to compare approximately the 
stabilities of the sections a to (a d,) and 
(a + d,) to (a + d) we must calculate the average 
velocity in the reverse flow section (a + d,) to 
(a +- d). 

From equations (2) and (3) 


d 
i : i da 
d 
d, 
d 
22, a | lia bx (d — 2)) 
did do). d 
d, 
$2, a (d,/d) (1 + d,/d) 


From equation (4) 


(2/9) 2, a 
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Dean's criterion may be applied to the reverse 
flow or outer section 


hes (d dy) (d — dy) 


From equations (9), (4) and (10) 


a 
297: 
22, a ‘ 4 AF d 


Assuming that TayLor’s criterion, equation 
(6), holds for the inner section a to (a + d,) 
and that the like an inner 
cylinder of radius a rotating in a stationary outer 
d,), the approximate 


(10) 
(11) 


section behaves 


cylinder of radius (a 
critical velocity for instability is 


Q.a = 41-1 
0 

dd 


Comparison of (11) and (12) shows that as 
2, increases instability will first occur in the 
d,). 


shown 


(12) 


section a to (a 
It will be that Tayor’s 

calculated for a full gap cannot be applied in this 

way for accurate prediction of critical velocity. 


criterion 


For partial reverse flow a similar prediction 
that instability will occur in the section a to 
(a + d,) can be made. 

In the comparison just made between the 
d,) and 
d) TayLor’s criterion was used 


stabilities of the two sections a to (a 
(a + d,) to (a 
assuming that there was a stationary solid 
surface at (a + d,); this is equivalent to a value 
of »« = 0 for cylinders of radii a and (a + dy) 
and would correspond to a value of S of the order 
of 3,400. It appears more reasonable to assume 
that the liquid flowing between (a + d)) and 
(a +- d) has an effect on the stability of the inner 
section and as a first approximation we could 
its effect is similar to that of the 


- d,) to (a + d) between 


assume that 
liquid in the section (a 
cylinders rotating in opposite directions at such 
a value of » that d,/d = 1/3. 

This would give a value of 2, characterized 
by S = 1,000 and the critical velocity would be 


Part I—Flow due to rotation of cvlinder 


1000 
reduced by 1000 or a factor of —— from that 
. 8400 1°85 


given by TayLor’s criterion for a stationary outer 
cylinder of radius (a + d,). 


It will be shown experimentally that the use 
of the value of S given by Fig. 9 for the particular 
value of d, /d considered gives good agreement and 


1S 


(13) 


Qa=1165 
dd 


when calculated from equation (8) 


S 1,000. 


using 


EXPERIMENTAL 


Calibration of apparatus 


Fic. 6. Diagram of experimental apparatus showing 


position of optical tank and end seals. 


The two cylinders of the apparatus (Fig. 6) were set 
horizontally. The inner cylinder was metallic and the 
outer made of Perspex. Built round the outer cylinder 
was an optical tank so that Section XX of the annulus 
could be observed. The speed of rotation of the inner 
cylinder could be varied over a wide range. Circumferential 
or end seals were used to limit the flow to a 25-cm length 
of cylinder. 

The apparatus was first calibrated by obtaining the 
critical speeds for liquids of different viscosities when 
the annulus was full. Indication of vortex formation 
was by dye injection when the liquid was water and by 
described by Hacerry [9] 
This effect 
observed with solution of viscosities as low as 8 centistokes 


self-indication 
with 


properties 
water/glycerine solutions. could be 
at room temperature. 

Values of 2, obtained for different values of viscosity 
are plotted in Fig. 7. 

By application of Tayvor’s criterion, equation (6), the 
effective value of d was found. 


d 0-867 cm, 


This value of d was used as the constant for the 


apparatus, 
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Fic. 7. Experimental critical angular velocity for Couette 
(full) flow and complete reversed flow. a 15-6 em, 
b 16-467 em (Calibrated by Taytor’s criterion). 


Completely Reversed Flow 


This experiment could be carried out quite simply 
with water, but with high viscosity glycerine solutions 
a thick film was carried over the top of the inner cylinder. 
Consequently, in order to have complete reversed flow 
with glycerine solutions a stationary barrier seal was 
inserted across the apparatus at YY, Fig. 6 

Values of critical speed were obtained for various 
values of viscosity and are given in Fig 


By 


the method of least squares 


Q. <= 35-275 
<= 1,244 


S was computed from equation (7a) using d)/d 1/3 


933 


The vortices produced were photographed (Fig. 8) 
using parallel transmitted light for illumination. Measure 
ment of the spacing of the vortices gave 


\d 5-45 em. 


This falls very nearly on the curve in Fig. 5 


Partially Reversed Flow 

The critical velocities were found for water for various 
free surface levels d’ between a and 6 and are plotted in 
Fig. 9. 


4 
Lt 
0-2 0-6 1-0 
ah 


Fic. 9. Experimental values of Q./v against reduced 
liquid depth d’/d for partially reversed flow. 
a 15-6 em b 16-467 om 
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Using the relation between d, d and d'/d 
given in equation (5) the curve in Fig. 9 was 
converted to one of 2. v against d, d. From this 


curve were calculated corresponding values of 


S from equation (7a). 

Fig. 10 gives Di Prima’s plot of Taytor’s 
results. On this the results of the present work 
are also shown. It is apparent that good agree- 
ment, in general, exists except at the value of 
d 


predicted values. 


d 1 where present results fall lower than 


7 
3400} + ; 'MNeoretica! valve 


3000 Refers to@=3-00 D= 4-035 
* Refers toad +55 C 4-035 
Refers to@ =3-80 O= 4035 
efers tod { 6-467 woter 
@ Refers tod = 15-6 5D =16-467 glycerine 


refers to a b 16-467 water 
@ refers to a 156 b 16-467 glycerine 


DISCUSSION 


Although a detailed mathematical analysis 
of the reversed and partially reversed flow 
problems would be desirable to establish the 
stability criterion, the simple laminar flow analysis 
together with the experiments performed have 
shown that the stability criterion may be obtained 
from the curve of S against (d,/d). 

Divergence of the author’s experimental values 
from those of TayLor as shown on Fig. 10, 
especially at d,/d = 1, are probably due to the 
inaccuracy of the apparatus used. The effective 
0-867 cm. was calculated from 
Taywor’s criterion using the experimentally 
determined value of Q./v for full flow. This 


value of d 


value of d when substituted in equation (7a) 
should have given a value of S = 8,400 if 
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Vortices produced by instability of completely reversed flow in 70 per cent glycerine solution. 


: 
| 


q 
> 
19057 
« 


The hydrodynamics of flow between horizontal concentric cylinders. Part I 


TayYLor’s criterion were exactly obeyed by the 
apparatus. In fact, as indicated in Fig. 10 for 
d,/d = 1, S = 2800, This error is probably due 
to the eccentricity of the cylinders, the inner 


being made from a 12in. i.d. copper pipe covered 


with 80 mesh wire gauze and the outer moulded 
in two halves from Perspex. No precautions 
were taken to reduce vibrations from the driving 


motor. 
CONCLUSIONS 


It has been established that the curve of S 
against d,/d provides the correct criteria for the 
instability of laminar flow between two concentric 
cylinders both when flow is of the Couette type 
and when reversed flow occurs. The calculation of 
d,/d using the approximate laminar flow equation 


is satisfactory. 
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NOTATION 
radius of inner cylinder 
radius of outer cylinder 
radial distance 
radius of lamina with zero velocity 
(b a) 
(To a) 
depth of liquid above inner cylinder 
(r a) 
angular velocity of inner cylinder 
angular velocity of outer cylinder 
29/2, = ratio of angular velocity of outer and 
inner cylinders 
tangential velocity 
kinematic viscosity 
pressure 
angle 
density 
critical angular velocity of inner cylinder 


REFERENCES 


Lams H. Hydrodynamics (6th Ed.). Cambridge (1932). 


Taytor G. IL. Phil. 


Trans. 


A 1922 223 289. 


Lewis J. W. Proc. Roy. Soc. A 1923 117 388. 


Synce J. 


Law C. C, 


Theory of Hydrodynamic Stability. 


L. Proc. Roy. Soc. A 1988 167 250. 


Cambridge (1955). 


Orrenneim A, K, and Hucues R. R. Industr. Engng. Chem. (Industr.) 1957 49 593. 
Dean W. R, Proc. Roy. Soc. A 1928 121 402. 
Di Pama R. C. Quart. Appl. Math. 1955 XIIT 55. 


Hacertry W. W. J. 


Appl. Mech. 1950 17 54. 


a 
b 
d 
do 
d 
$2) 
(Jo 
VOL. 
957/58 
[6] 
[7] 
(8) 
221 


Pergamon Lid 


Chemical Engineering Science, 1958, Vol. 7, pp. 222 to 228. 


On the mechanism of interfacial resistance to mass transfer 
liquid-liquid extraction 


G. Boye-Curistensen and G. 


The Technical University of Norway, Chemical Engineering Laboratory, Trondheim, Norway 


(Received 16 May 1957; revised 30 August 1957) 


Abstract—The additional resistance to the transfer of o-nitrophenol from an aqueous solution 
to falling drops of carbon tetrachloride caused by the addition of the surface-active agent sodium 
oleyl-p-anisidine sulphonate has been determined, and the results compared with those previously 
reported with iodine as the diffusing solute. 

The results are adequately described by the adsorption isotherm previously developed, 
but the adsorption equilibrium appears to be more closely approached with smaller drops than 
with large ones. The equation of state for the adsorbed molecules, however, applies irrespective 
of the degree of approach to equilibrium, 

The percentage reduction in the rate of mass transfer is shown to be 66-69 per cent for this 
system both with iodine and o-nitrophenol as diffusing solutes and for drop sizes between 7-5 
and 18-5 mm®. It is concluded that there is no specific chemical or physical interaction between 
the adsorbed surface-active molecules and the diffusing solute, and that the results favour a 
hydrodynamic type of mechanism for interfacial resistance. 


Résumé—Dans le cadre de l'étude de l'échange entre l'o-nitrophenol en solution aqueuse et une 
pluie de gouttes de tétrachlorure de carbone, les auteurs ont étudié la résistance supplémentaire 
causée par l'adjonction de l'agent tensio-actif oleyl-p-anisidinesulphonate. De plus, les valeurs 
obtenues ont été comparées avec celles déja obtenues en employant l’oide comme soluté diffusant. 


Les résultats sont correctement rapportés par lisotherme d’adsorption déja deécrit ; 
toutefois l'équilibre parait étre plus rapproché avec des petites gouttes qu’avec des grandes. 

L'équation d'état pour les molécules adsorbées s'applique cependant, sans tenir compte du 
degré d’approche de l’équilibre. 


La réduction affectant l'échange massique, exprimée en pour cent, est de 66 — 69 pour les 
deux solutés diffusant utilisés—liode et l'o-nitrophénol—et pour des gouttes de 7,5 & 18,5 mm3. 
Les auteurs en concluant qu'il n'y a pas d'action chimique ou physique spécifique, entre les 
molécules actives de la surface d’adsorption et du soluté diffusant. Ces résultats les aménent a 


conclure en faveur d'un mécanisme du type hydrodynamique pour la résistance interfaciale. 


Zusammenfassung—Es wurde der zusitzliche Widerstand auf die Cbertragung von o-Nitro- 
phenol aus einer wissrigen Lisung auf fallende Tropfen von Tetrachlorkohlenstoff bestimmt, der 
durch die Zugabe des oberflichenaktiven Stoffes Natrim-Oleyl-p-Anisidinsulfonat hervorgerufen 
wurde. Die Ergebnisse wurden mit den friher mitgeteilten verglichen, bei denen Jod der diffun- 
dierende Stoff war. 


Die Ergebnisse werden hinreichend durch die friiher mitgeteilte Adsorptions-isotherme 
beschrieben, jedoch scheint das Adsorptionsgleichgewicht besser bei kleineren als bei grésseren 
Tropfen erreicht zu werden. Die Zustandsgleichung der adsorbierenden Molekiile scheint 
unabhangig von der Anndherung an das Gleichgewicht zu sein. 


Die prozentuelle Verringerung der Grésse der Stoffiibertragung ergibt sich zu 66 bis 69% 
fiir dieses System mit Jod und ¢-Nitrophenol als diffundierende Substanz und fiir Tropfengréssen 
zwischen 7,5 und 18,5 mm*, Daraus kann man schliessen, dass es keine chemische oder physikal- 
ische Wechselwirkung zwischen den adsorbierten oberflichenaktiven Molekiilen und der diffun- 
dierenden Substanz gibt und dass die Ergebnisse auf einen hydrodynamischen Mechanismus fiir 
den Grenzflichenwiderstand hindeuten. 
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INTRODUCTION 
Tue nature of the additional resistance to mass 
transfer in liquid-liquid extraction caused by 
the 
been discussed by several authors. 


substances has 
For 
venience, the term interfacial resistance is used 
for this phenomenon. The term implies that the 


presence of surface-active 


con- 


resistance is caused by surface-active molecules 
not be 
taken to give any indication of the responsible 


adsorbed at the interface, but it must 


mechanism, The explanations advanced can be 
classified either as barrier theories or as hydrodyna- 
mie theories. The former category is based on 
the assumption that surface-active molecules 
adsorbed at a liquid-liquid interface form a kind 
of barrier which interferes with the passage of 
diffusing solute molecules, whereas the latter 
presupposes that the adsorbed molecules are 
somehow capable of modifying the hydrodynamic 
conditions in the vicinity of the interface. 
Hurcurinson [1] found that the resistance was 
highly specific. He favoured the barrier theory, 
and suggested that the action is not mechanical, 
but involves forces of a physical or chemical nature. 
West, Beck 
McGrecor [2] found that the rate of extraction 
of acetic acid from benzene by drops of water 


Ropinson, MorGeNTHALER, and 


was greatly reduced by traces of plasticizer 
extracted from Tygon tubing in the benzene feed 
line. They concluded that this effect is associated 
with a reduction of the internal circulation in the 
aqueous drops. In a later paper West, Herrman, 
Cuonc and Tuomas [3] studying the effect of 
impurities in the the 
highly tentative suggestion that the surface- 
active substances act as a kind of interfacial 


benzene phase made 


barrier. 

Garner and Hare [4] adopted the view of 
Hutcuinson that the interfacial resistance is 
caused by physical or chemical interaction 
between the diffusing solute and the surface- 
active molecules at the interface. GARNER and 
Hae [4] and Garner and SKELLAND [5] have 
shown that additions of surface-active agents 
reduce the rate of internal circulation in the 
drops, but that diffusing solute can counteract 
this effect to some extent. They found, however, 
that in extracting diethylamine from toluene by 
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drops of water no internal circulation could be 
observed even in the absence of surface-active 
agent. On this basis they rejected the hydrodyna- 
mic theory in favour of a barrier type of theory. 

According to Lewis [6] transfer of turbulence 
across an interface constitutes an important 
element in mass transfer, and he explains inter- 
facial resistance as a result of a reduction in the 
transfer of turbulence 
surface-active molecules. 


caused by adsorbed 
The effect of adsorbed 
molecules on the spontaneous interfacial turbu- 
lence described by Lewis and Prarr [7] must 
also be included among the possible explanations. 
The situation has recently been discussed by 
LinpLanp and Tersesen [8], but the experi- 
mental material so far available has been in- 
sufficient for any firm conclusions to be drawn. 
Some advances have, however, been made along 
more general lines. GARNER and Hate [4] have 
that 
interfacial 
and Me.uus and Tersesen [9] have concluded 


shown there is no parallelism between 


resistance and _ interfacial tension, 
that interfacial resistance and interfacial tension 
are two different manifestations of adsorption, 
unrelated to each other except for the fact that 
Whereas 
Garner and Hace [4] did not find any direct 


they both originate in adsorption. 


connection between adsorption and_ interfacial 
resistance, the latter authors developed a theory 
of interfacial resistance, assuming, as a first 
approximation, proportionality between the two 
phenomena. 

Extensive determinations of interfacial resis- 
tance are available for the system water - carbon 


tetrachloride with iodine as diffusing solute 
and sodium oley!l-p-anisidine sulphonate as surface 
active agent [8]. 
particularly suitable for studies of interfacial 
resistance, and it is the purpose of the present 
investigation to extend this work by substituting 
o-nitrophenol for 
Dissolved molecules of iodine and o-nitrophenol 
have very different physical and 
properties, and any physical or chemical forces 
operating at the interface must be expected to 
interfere differently with the passage of the two 


This system appears to be 


iodine as diffusing solute. 


chemical 


substances. 
with 


A comparison of results obtained 


these substances should therefore yield 
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additional information about the nature of 


interfacial resistance. 


EXPERIMENTAI 

The apparatus and experimental procedure has 
been described in detail elsewhere [8]. Some 
difficulties were experienced in maintaining 
constant drop volume throughout a run, but 
capillaries drawn from thin-wall glass tubing 
gave satisfactory performance. An observation 
of GARNER and SKELLAND [13] that the time of 
drop formation is important was confirmed, 
and their suggestion to work with fairly low 
drop frequencies was adopted. The temperature 
was 25°C as before. and the concentration of 
o-nitrophenol about 0-020 ¢ 100 ml aqueous 


phase. 


MATERIALS 


The carbon tetrachloride was of technical 
quality, and fractionated in a column of 5-10 
theoretical plates with a reflux ratio of about 6: 1. 
The fraction used distilled between 76 and 77°C 
and had a refractive index of n> 14568 as 
compared with the value of 1-4571 given in the 
International Critical Tables. When a sample 
was distilled in an analytical column of 60-80 
theoretical plates at a reflux ratio of 99:1 at the 
beginning and end of distillation and 49:1 in 
the middle, about 97 per cent went over at 
76-8°C (760 mm Hg). The forerunnings and the 
residue had boiling points very close to this value 
and approximately the correct refractive indices, 
showing that the solvent phase consisted of 
practically pure carbon tetrachloride. 

The o-nitrophenol was Merck Pro Analysi. 
It is an acid with a dissociation constant of about 
7-5 x 10° and in order to reduce dissociation to 
a negligible value the pH of all aqueous solutions 
was adjusted to 3-15 by making the solution 
0-002 N with regard to citric acid (Merck Pro 
Analyst). 

The sodium oleyl-p-anisidine sulphonate was 
purified and stock solutions prepared with sodium 
sulphate and sodium chloride as described 
previously [10]. Stock solutions older than one 


week were discarded. 
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Mretuops or ANALYSIS 


Solutions of o-nitrophenol in water at 
pH = 3-15 and in carbon tetrachloride were 
analysed photometrically using a Beckman DU 
quartz spectrophotometer. The readings were 
taken at the absorption maximum which was 
found to be located at 3500 A for both aqueous 


and solvent solutions. 


ResuLTs AND Discussion 

Experiments were carried out with drops 
varying from 5 to 22 mm" in volume and with 
bulk concentrations of surface active-agent up to 
3-5 10+ ¢g ‘100 ml of aqueous phase. Larger 
drops gave poor reproducibility, and at higher 
concentrations of surface-active agent the drops 
coalesced rather slowly in the receiver at the 
bottom of the column. 

During the earlier part of the work it was not 
possible to maintain constant drop volume from 
one run to another, and it was necessary to 
determine the percentage approach to equilibrium 
p as a function of drop volume for each concentra- 
tion of surface active-agent. The volumes were, 
however, kept as close to 7-5 and 18-5 mm® as 
possible, The values of p for these two selected 
volumes were then determined by graphical 
interpolation. This applies also to the rate of 
fall of the drops. The results are given in Table 1 
together with the corresponding total and inter- 
facial resistances R, and Ry, calculated on a time 
and extract or solvent basis as described previously 
[8]. In addition, the resistances are given also 


on a raflinate or aqueous basis : 
Ry = and Ra, = Ryy/m 


The distribution coetflicient for o-nitrophenol 
between water and carbon tetrachloride was 
required for the calculations, It was determined 
by alternately shaking and leaving in a thermostat 
at 25°C a glass-stoppered bottle containing 50 ml 
of 0-002 N aqueous citric acid and 50 ml of 
carbon tetrachloride to which was added a 
suitable amount of o-nitrophenol. Within the 
range of 0-011 to 0-022 g of o-nitrophenol per 
100 ml of aqueous phase the result was : 
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sulphonate g/100 mL aqueous phase 


Interfacial resistance coefficients on a raflinate or aqueous basis with o-nitrophenol as diffusing solute. 


O Drop volume 7-5 


m = 126-5 + 2-5 
mg of o-nitrophenol per ml of solvent phase 


mg of o-nitrophenol per ml of aqueous phase 


The distribution equilibrium is so greatly in 
favour of the solvent phase that the resistance to 
mass transfer can be considered to be located 
almost entirely on the aqueous side, as was the 
case also with iodine as diffusing solute. 

The interfacial resistances have been plotted 
on Fig. 1 where the curves represent the best 
fitting equation of the type developed previously 
[8] : 

i 
R ol (1) 
ix 
Values of the constant K in equation (1) are 
given in Table 2 together with the limiting 


interfacial and total resistances RR; and R,, 


based on present work and previous results with 


iodine as diffusing solute [8]. The values have 
been calculated on a raflinate or aqueous basis as 
well as on an extract or solvent basis. The 
limiting values of the resistances represent 
the average for the concentration range 
Cy = 10 — 85 x 10° ¢/100 ml. The constants 
in the corresponding isotherm : 


@ Drop volume 18-5 mm? 


kK 


Ss isotherm (1 


are included also. 


6 7 8 
Concentration of sodiurn oley+p-anisidine 
sulphonate g/1OO ml aqueous phase 


Fic. 2. Adsorption isotherms. 
o-nitrophenol, drop volume 7-5 mm® 
o-nitrophenol, drop volume 18-5 mm® 
Iodine, average for drop volumes 4-2, 9-5 and 
22-3 mm® 


It can be seen from Fig. 1 that the theoretical 
equation is capable of describing the results 
within the limits of experimental error. The 
results, which are more accurate than the 
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Table 1. Summary of results 


7 


Cs x 108 Ry x 10 Rg, = 16% | Rp 
(2/100 mil) (cm /sec) m Cra (sec /em) (sec /em) (sec/com) (sec 
Drop volume: 75mm, equivalent radius: 1-21 mm 
0 18-8 | 8-12 5-9 0 ‘7 0 
170 1-84 112 5-3 42 
2 16-8 1-46 14-4 8-5 114 67 
3 16-6 1-27 16-7 10-8 132 85 
5 16-5 1-23 17-3 137 
10 16-5 1-23 17+ | 11-4 137 90 
25 16-3 1-23 17-6 | 11-7 139 | v2 
35 16-3 1-25 ll-4 137 90 
Drop volume: 185mm, equivalent radius: 1-64 mm 
0 22-0 2-03 | 58 0 46 0 
1 19-7 1-34 9-8 | 40 78 | 32 
2 19-3 1-04 13-0 72 103 57 


2 0-87 15-7 9-9 
4 2 17-6 11-8 
5 19-2 0-76 18-0 12-2 145 
1v-2 0-73 18-7 12-9 148 102 
0 O75 18-2 12-4 07 
. 19-2 0-73 18-7 | 12-9 148 102 
1 18-8 13-0 


12-8 


18-6 


previous ones, also show the effect of drop size active agent the fraction of the interface covered 
This 


becomes particularly clear from an inspection of 


by adsorbed molecules is considerably larger for 


on the action of the surface active-agent. 


drop volumes of 75mm" than of 18-5 mm’, 
With the smaller drops the interface is practically 


the slope of the adsorption isotherms in Fig. 2. 


bulk 


concentration of the surface- completely covered when the bulk concentration 


For a given 


Limiting resistances and values of the constants in the equations (1) and (2) 


Table 2. 


Reference Present work Lindland and Terjesen 


Diffusing solute o-nitrophenol Iodine 


Drop volume (mm?) 18-5 75 22-3 0-5 4 


seesem 18600 17400 11900* 10800 
Rp.w, sec/em 12800 11400 Average value 7400 
sec/em 147 138 130 120 
see/em 101 91 Average value 82 
7 Kg, gem/see 100 ml 26 x 10°° 21 x 10° Average value 5-4 x 10°° 
. Kp, gem/sec 100 ml 33 x 107 26 x 107 Average value 4-9 x 10-7 
e Kigotherm: £/100 ml 3-3 x 10°5 24 x 10° Average value 40 x 10°5 
| 


*In ref. [8], in the second column on page 8, these values are incorrectly called the overall resistances without 


addition of surface-active agent. 
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is about 45 « 10° ¢/100 ml of aqueous phase 
whereas about 5-5 x 10° g/100 ml is necessary 
with the larger drops. The adsorption equilibrium 
is thus more closely approached the smaller the 
drops. 

The exponent in the denominator of equations 
(1) and (2) is an expression of the sharpness of 
the transition between the rising and _ the 
horizontal branch of the curves. 
from Fig. 1 that the value 0-1 fits very well in 
both cases. It has previously been shown [9] 
that the exponent, but not the constant, in the 
isotherm enters into the equation of state for the 
adsorbed film and therefore must be expected 
to be invariant with respect to the degree of 


It can be seen 


approach to adsorption equilibrium. 

The average isotherm obtained previously with 
iodine as diffusing solute exhibits a somewhat 
slower approach to equilibrium, It must, however, 
be remembered that the surface-active agent in 
that case had added two molecules of iodine to 
the double bond, and therefore really it was a 
slightly different substance. 

For the purpose of comparing the different 
results the interfacial resistance R,; has been 
expressed as per cent of the total limiting resis- 
tance R, and plotted in Fig. 3, The total 
limiting resistance given in Table 2 has been 
used as a basis rather than the resistance without 
addition of surface-active agent because the 
former appears to be more accurate, particularly 
in the case of the previous work. The resistance 
without addition of surface-active agent is rather 
sensitive to traces of impurities. It can be seen 
from Fig. 8 that for concentrations of surface- 
active agent below 6 x 10° g/100 ml there is a 
small, but systematic deviation from the fully 
drawn curve representing the average of all 
observations. 

For higher concentrations, on the other hand, 
the interfacial resistance amounts to 66-69 per 
cent of the total limiting resistance, irrespective 
of diffusing solute and drop size. This is equiva- 
lent to stating that the percentage reduction 
in the rate of mass transfer caused by additions 
of surface-active agent greater than 6 x 10° 
g/100 ml is 66-69 per cent. Since the two 


solutes differ considerably in physical and 
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8 
x 
3 
> 
2 4 6 0 12 14 6 18 20 22 24 26 28 
Concentration of sodium oleyl-p-anisidne sulphonate 9/100mt 
aqueous phase 
Fic. 3. Interfacial resistance coefficients R; expressed 
as per cent of the total limiting resistance coefficients R,. 
VY o-nitrophenol, drop volume 7-5 mm? 
A e-nitrophenol, drop volume 18-5 
0 Iodine, drop volume 4-2 mm® 
O Iodine, drop volume 9-5 mm? 
@ Iodine, drop volume 22-3 mm’ 
chemical properties the retardation of mass 


transfer therefore appears to be caused by a 
general mechanism not involving any specific 
chemical or physical interaction between the 
diffusing solute and the adsorbed surface-active 
molecules. 

The observation that a change from one solute 
to another affects the interfacial resistance and 
the normal resistance in approximately the same 
proportions can be explained most easily by 
assuming that interfacial resistance is the result 
of modifications to the normal mechanism of mass 
transfer rather than an introduction of a new and 
additional mechanism. Such a modification to 
the normal mechanism of mass transfer can take 
place through changes in the hydrodynamic 
conditions in the vicinity of the interface caused 
by the adsorbed surface-active molecules. 

The limiting total resistances Rp, are fairly 
accurate for iodine as well as for o-nitrophenol. 
It is of interest to consider the ratio of these 
resistances for the two different diffusing solutes, 
but for approximately the same drop size. The 
result of such a comparison is given in Table 3. 

It is generally agreed that liquid phase diffusivi- 
ties are inversely proportional to the molar volume 
of the solute to the power of 0-6. The molar 
volumes of iodine and o-nitrophenol have been 
calculated according to the method of WILKE [11] 
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Table 8. Comparison of limiting resistances with 


Drop size, mm® : 


For iodine 22:3 9-5 
For o-nitrophenol 18-5 75 

; 
lodine_ 0-88 0-87 


nitrophenol 


and Witke and Cuane [12], and found to be 
71-5 and 130-7 em*g mole respectively. The ratio 


of the diffusivity of o-nitrophenol to that of 


iodine is therefore 0-70. On the assumption that 
the difference between the experiments with the 
two solutes is represented entirely by the two 
diffusivities the results can be expressed by the 


equation : 0-38 (3) 


CONCLUSIONS 


(1) Both the additional resistance caused by 
the presence of the surface-active agent and the 
normal resistance are affected in the same pro- 
portions when o-nitrophenol is substituted for 
iodine as diffusing solute. This indicates that the 
additional resistance is a result of a modification 
to the normal mechanism of mass transfer rather 
The 


results thus favour a hydrodynamic type of 


than an introduction of a new mechanism. 


theory. 

(2). The rates of mass transfer is proportional 
to the diffusivity to the power of 0-38. 

(3) The absorption equilibrium is more closely 
approached with smaller drops than with large 


ones. 


Hurcuinson E. J. Phys. Chem. 1948 52 897. 


(Industr.) 1951 43 234. 
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NOTATION 


Ay effective cross-sectional area of 


adsorbed surface-active molecules em?/g mol 

ag = effective cross-sectional area of ad- 
sorbed surface-active molecules A? ‘molecule 

Cg = concentration of o-nitrophenol in 
solvent phase mg /em3 

equilibrium concentration in sol- 

vent phase mg /em8 

Cr = concentration of o-nitrophenol in 

8 


aqueous phase mg/em 


average of initial and final con- 
centration in aqueous phase mg/em3 
Cg = concentration of  surface-active 
agent in aqueous phase 
D = diffusivity 


m distribution coeflicient : 


z/100 ml 


em? /sec 


mg o-nitrophenol per ml solvent phase 


mg o-nitrophenol per ml aqueous phase 
Rg total resistance to mass transfer on 
a time and extract or solvent basis sec 


c/em 
Ry; = interfacial resistance on a time and 
solvent basis see ‘cm 
Rp = total resistance on a time and 
raflinate or aqueous basis sec /em 
Rp; = interfacial resistance on a time and 
aqueous basis sec /cm 
Ry limiting resistance i.e. resistance 
reached when further additions of 
surface-activeagent caused no fur 
ther increase sec /cm 
v = drop velocity em /sec 


I’ = surface concentration of surface- 
active agent 


fom? 
g mol/em 
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Abstract 
and energy for a rapid reaction occurring in the laminar boundary layer on a flat plate. The 


A simple, exact (boundary-layer) solution is obtained to the equations of diffusion and 


velocity distribution is assumed to be given by the Blasius solution to the Prandtl boundary- 
layer equations. Expressions are derived for rates of interphase transfer, and concentration and 


temperature distributions. Applications are to problems of absorption, extraction, dissolution 


and combustion. 


Résumé—C’est la solution exacte et simple des équations de diffusion et d’énergie qui est ici 
obtenue pour une réaction rapide ayant lieu dans la couche limite laminaire d'une plaque plane. 
Il est admis que la distribution des vitesses est celle fournie par la solution de Blasius aux 
équations de la couche limite de Prandtl. Ceci permet de déduire des expressions pour les taux 
Des 
applications peuvent étre faites aux problémes d’absorption, d’extraction, de dissolution et de 


déchanges dans linterphase, et la répartition des températures et des concentrations. 
combustion. 
Zusammenfassung—Fiir die Gleichungen der Diffusion und der Energie bei einer schnellen 


Reaktion in Platte wird exakte 
(Grenzschicht—) Lésung angegeben. Die Geschwindigkeitsvertcilung wird gemiiss der Blasius- 


der laminaren Grenzschicht einer ebenen eine einfache, 


Fiir den Austausch zwischen 
den Phasen und die Verteilungen von Konzentration und Druck werden Ausdriicke abgeleitet. 


Lésung der Grenzschichtgleichungen von Prandtl angenommen. 


Anwendungen stellen dar: Probleme der Absorption, Extraktion, Auflésung und Verbrennung. 


INTRODUCTION 
the last 
amount of information, both theoretical and ex- 


to heterogeneous catalysis in a flow system. 
KISER Hoe.scuer [7] describe a water 
tunnel constructed to facilitate hydrodynamically 


DURING two decades a considerable and 


perimental, has become available on the dynamics 
of real fluids. Much of this information is sum- 
marized in the monographs on turbulent flow by 
Barcuetor [1] and Townsenp [15] and on 
[10]. 
Recently, investigators have sought to apply 
these results to moving fluids in which chemical 


boundary-layer theory by ScHiicutTinc 


reactions are taking place, systems of great 
practical importance. Thus Corrstn [8] has begun 
the development of a statistical theory for chemical 
reaction in a turbulent fluid, and Cuampre and 
Acrivos [2] have applied boundary-layer theory 


controlled studies of catalytic surface reactions. 

Another important class of chemical changes in 
flow systems is the rapid, irreversible second- 
order reaction. Examples are the absorption of 
carbon dioxide by dilute solutions of base [4] 
or the dissolution of solid acids in basic solutions 
[5, 6]. A similar model is used in combustion 
calculations for the burning of single fuel drops 
[8, 14] but here the picture is further complicated 
by large heat effects. 

As pointed out by SuHerwoop and Picrorp 
[13], two approaches are commonly employed in 
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making quantative predictions for these systems. 
In one, a modification of the old Nerst-Brunner 
theory, it is postulated that a stagnant film 
exists next to the interface. The component which 
is dissolving, evaporating, or being absorbed 
diffuses across the phase boundary to a point 
within the film where it meets the other reactant 
diffusing in from the main body of the liquid or 
gas. Reaction is instantaneous. In the steady 
state, the reaction zone takes up a position which 
depends on the diffusivities and concentrations 
of the components. 

Alternatively, absorption in packed columns 
has been pictured as an unsteady state pheno- 
menon in which packets of liquid contact the 
interface for a short time and then return to the 
main bedy. During the period of contact, un- 
steady molecular diffusion to the reaction zone 
(as defined above) takes place but the position of 
the reaction zone changes with time. 

These models, while easily pictured and 
mathematically tractable, are not faithful re- 
presentations of any real flow system. As a guide 
to the influence of convective effects it would be 
useful to have theoretical predictions for an 
actual system, well defined hydrodynamically. 
It is shown below that with some not too res- 
trictive assumptions it is possible to obtain a 
simple, exact (boundary-layer) solution to the 
problem of rapid chemical reaction in the laminar 
boundary layer on a flat plate. 


ANALYSIS 


As an example, the steady state dissolution of a 
slightly soluble solid, such as benzoic acid, in a 
flowing dilute solution of a base may be con- 
sidered. The geometry is that of a flat plate of 
acid at zero incidence to the flow. The process is 
practically isothermal, and the properties of the 
fluid are independent of position. The rate at 
which dissolution takes place as a function of 
distance from the leading edge of the plate and 
concentration of reactant in the mainstream (at 
infinity) will now be discussed. 

For low rates of dissolution, the normal 
velocity at the interface is small. Provided the 
Reynolds number is_ sufficiently large, the 
velocity distribution is given by the Blasius 
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solution of the Prandtl boundary-layer equations 
ff" + 2f" =0 (1) 


where f’ (») represents the (dimensionless) 
boundary-layer velocity and = U/va. 
Values of f and its derivatives are tabulated in 
Scuuicutine [11)}. 

Substance A dissolves from the plate and 
diffuses to the reaction zone while B diffuses into 
the reaction zone from the main body of the 
solution. In the reaction region, which is con- 
sidered to be of negligible thickness, the two 
substances meet, and the rapid, irreversible 
reaction : 

aA + bB - products (2) 
takes place. In the region between the reaction 
zone and the surface the diffusional process is 
described by the expression :— 

oy oY 

Here the term representing diffusion in the x 
direction is neglected in accordance with the 
usual boundary-layer assumptions. Similarly in 
the region outside the reaction zone :— 

oc 

oY oY 

Transforming to the » co-ordinate system, it is 
easily shown [12] that these equations become :— 


dts , (5) 
S4 d7* 2 dy 
and 
Sp 2 dy 
where 
= 
C4= Cg = (7) 
Cus Crm 


The boundary conditions imposed on the system 
are as follows :— 


7 = "pes = Cy =O (8) 


Stated in words, it is assumed that the reaction 
takes place along a surface of constant 7 = n,z. 


= 
1957 
At = 0, ¢,=1 
= Cr = : 
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This assumption is necessary to preserve the 
similarity postulated in the transformation to the 
7 co-ordinate system. The final condition necessary 
to fix the position of y, results from the assump- 
tion that the materials enter the reaction zone in 
stoichiometric proportions : 


At DR 
(9) 


(9a) 


Since f = — 2f'"/f’’, equations (5) and (6) can 
be integrated between limits to give :— 


. 


dy 


0 (10) 


f° dy 


” 
. 


dy 


"R 


| dy 


R 


Substituting in (10) for ¢4 and cy, and re- 
arranging the following equation is obtained : 


S4) \0-332, 1 — Oar 


using the approximation (12) :— 


bC 4s 


(0-332) 


0-332 S83 (18) 


The parameter @, is defined by the expression: — 


(n, S) = 


and is plotted by Scuitcurinc [12] after 
PounLHausen [9]. From equation (10) the local 
Sherwood number is : 


N, x 


0-332 S 
dy 


Ov AR 


and the mean Sherwood number for the entire 
plate up to 2 is given by :— 


l 
— [0-664.8 16) 


The quantity 1/1 — 0,4, represents a correc 
tion factor for reaction when applied to the term 
0-664 S ,* R,}, the Sherwood number for pur 
Since 0 <1, 
rate for simultaneous reaction will always be 


diffusion. the mass transfer 
greater than that for diffusion alone. 

For a given system, with S,, S, and C 4g fixed, 
the reaction factor can be obtained as a function 
of main stream concentration in the following 
way: A value of nz is assumed and @, 4p, 
and f, read off the appropriate curves and 
tables [11, 12]. The 
calculated immediately and the corresponding 


reaction factor can be 


value of aCyy/bCys can be obtained from 
Fig. 1 the 
a calculation of this type taking S 


results of 
= 900 for 


equation (12). shows 
benzoic acid and 524 for sodium hydroxide in 
water at 25°C. The values of S are based on the 
data given by Kino and Bropte [6]. 
For the 
= and from equation (12) :— 


important case of S, = Sp, 4p 


l 
al py 
bCas 


and 


(18) 


It is interesting to note that the reaction factor 
that 
obtained for the simplified model of constant 
film thickness without (13). The 
factor for S., = Sy is also shown in the Figure. 

When Cyy,/C 4s is very large the reaction zone 
approaches the surface, fp = 0-332 and 
1. Equation (12) can be rearranged to 
give 


1 +{(aCyy)/(bC4s)} is identical with 


convection 
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| | dy 
(14) 
| dy | 
0 
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: was found to be negligible. 
Pr 
fad, For a gaseous reaction S ~ P and equation (24) 
1 >» > 0 (21) becomes : 
; 
R 
| AH 2D 4s/k = (Tp — Ts) 
— Ty) (25) 
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ac py (19) The position of the reaction surface np can be 
Qian Cas \Sp 


calculated from the flow-diffusion equation in a 


manner described in the previous section. Its 
rhis explains the linearity of the curves in Fig. . , 
: value depends on the temperature distribution 
for large values of C py, /C 4s. 

- through the solubility or vapour pressure relation- 


—-- - ship Cy = F (Ts). 
| The temperature of the reaction zone is obtained 
| by assuming that the heat generated at » is 
by é 4 conducted away from both sides of the zone 
o 4 Ze af according to the expression : 
4 on 


Fic. 1. Reaction factor as a function of concentration ratio, AH,D 


EXTENSION TO NoN-IsoTHERMAL SYSTEMS (Te — Te) tan + (Te T tye (28a) 


S “AR 


This approach can be extended to the steady ; 
state non-isothermal case with constant fluid (Since the reaction zone ts of vanishing thickness, 


properties to permit calculation of the tempera the heat convected in and out is equal), Re- 


ture of the reaction zone and the distribution writing equation (2 sa) in terms of ¢ Pe 
near the surface. The energy equivalent to 


equation (5) Is : 


1 {dT 
P dy 2 dy 


MyD Cas 
kil Ov ar) P 


SA 
fro 


0 (20) Tr—-Ts . Tr 


assuming that dissipation by viscous friction is 

T, and where = 9, XR) and Oy = 9, (P, np). 

the If », is known from the diffusion calculation 


that the surface and mainstream temperatures 


negligible. If the surface temperature 


main stream temperature Ty) are constant, 


temperature 7', of the reaction zone | t also be can be obtained from equation (24). Assuming 


constant to preserve similarity. Under thes 


conditions, the solution to equation (20) for each were both 25°C, a caleulation of this type was 


side of the reaction zone is identical with the form ™#de for the system benzoic acid-sodium 


for the concentration distributions : hydroxide solution. For dilute solutions of the 


base, the temperature rise of the reaction zone 


Flame front temperatures for combustion 

- processes can be estimated from equation (25) 
uP 

| f° B dy provided that mean fluid properties can be used 


in the equations of motion, diffusion and energy. 
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APPLICATION TO OTHER GEOMETRIES 


The 
laminar boundary layer on a flat plate can be 
extended to for which 


similar (one space parameter) solutions to the 


* exact ’’ method developed above for the 


other laminar flows 


boundary-layer equations exist. Examples are 
flow past a wedge and free jet flow. Application to 
turbulent S\ stems (such as free jets) may also be 
possible provided the reaction takes place over a 
sufliciently narrow region. 

For shapes such as eylind rs or spheres which 
the 


do not permit similarity transformations, 


von Karman integral method can be applied. 
Here it is necessary to assume the existence of 


two boundary layers; at the edge of the inner 


layer, the boundary conditions are given by 


equations of the type (8) and (9) and a polynomial 
form can the in 
each layer. A calculation of this type was made 


for the flat plate and the results were found to 


be assumed for concentration 


compare well with those of the exact solution. 
Over a large Reynolds number range most of 
the transfer for such shapes as cylinders and 
laminar boundary layer 
the 


Sherwe 


spheres occurs in a 
back 
point ; the 


diffusion is given by 


Ny 


extending from forward stagnation 


number for pure 


here 


The quantity A is a constant which depends on 
the geometry but is usually near 0-664, the flat 
plate value. Thus, as a first approximation, tt 
that the 


factor for other geometries will be of a form 


seems reasonable to assume reaction 


similar to that for the flat plate. 
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NOTATION 
stoicheiometric coeflicients of chemical reaction 
boundary-layer concentration parameters for 
reactants 4 and B, dimensionless 
boundary-layer concentrations of reactants 4 
and B, moles/1 
surface concentration of A, moles/1! 
mainstream concentration of B, moles 
heat of fluid, cal/g ¢ 


of 


ifi 


diffusivities reactants A and B in fluid 


cn.~, sec 
dimensi mnless 


al 
fluid 


boundary-layer velocity 
standard heat of reaction, k vy mok 


thermal conduc tivityal eal /see (em* 


mies 
Les 
Prandtl! number for fluid, cp p ‘v/k 
distance Reynolds number, Ua/: 

Dg 


Schmidt nu mber for reactant B. Dp 


sce 


il mass transfer coeflicient, cm 
in Sherwood number, 2Ny dimensionless 


il Sherwood number, k,.7/ D4, dimensionless 


S-hmidt number for reactant 


boundary-layer temperature parameters, di- 


mensionless 

temperature of surface, ¢ 

temperature of reaction zone, ( 

temperature of fluid main stream, “¢ 
boundary-layer velocity in # and y directions 
respectively, em /se¢ 

main stream fluid velocity, em /sec 
direction parallel to plate 
direction perpendi« ular to plate 
sec 


kinematic viscosity, 


dimensionless coordinate, y 4/ U /av 
density of fluid, g /« m® 


denote differentiation with 


‘ R 


Notes :—Primed quantities 
respect to n: 


at » 


the subscript refers to values taken 


oR 
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Countercurrent gas-liquid flow 


Stability of a spherical liquid film in an infinite airstream 


J. E. Howxins and J. F. Davipson 


Chemical Engineering Department, Cambridge 


( Received 8 August 1957) 


Abstract—A liquid film runs symmetrically under gravity over a string of spheres suspended in 


an upward stream of air. At a certain air velocity a standing wave forms on the film covering the 


lowest hemisphere. The wave is formed at about 30° to the vertical, where the pressure gradients 


within the airstream, due to its acceleration over the sphere, are a maximum 


Experiments with spheres of 1-0 and 1-5 in. diam. suspended in a 5-75 in. tube gave good 


agreement with theory. Wave formation of this kind is thought to be responsible for loading in 


packed towers of high voidage 


Résumeé-—Les auteurs considérent un film liquide s’écoulant symétriquement par gravité sur 


un chapelet de sphéres suspendues dans un courant d’air ascendant. Une onde stable se forme 


sur le film couvrant 'hémisphére inférieure, pour une certaine vitesse de l'air. Cette onde est placée 


& 30° environ sur la verticale 4 l'endroit of les variations de pression dans le courant d’air, ddes A 


l'accélération sur la sphére, sont maximum. 


Des sphéres de 1 A 1,5 pouce de diamétre suspendues dans un tube de 5,75 pouces donnent 


des résultats en accord avec la théorie, C'est sans doute cette onde que détermine le point de 


charge dans les colonnes 4 garnissage ayant un grand degré de vide. 


Zusammenfassung-—Ein flissiger Film liuft symmetrisch unter dem Einfluss der Schwerkraft 
iiber ein Band aus Kugeln, das in aufwirts stroémender Luft aufgehingt ist. Bei einer bestimmten 
Geschwindigkeit bildet sich auf dem Film, der die unterste Halbkugel bedeckt, eine stehende 
Welle. Die Welle bildet sich bei cinem Winkel von etwa 30° von der Senkrechten, wo die 
Druckgradienten des Luftstromes, hervorgerufen durch seine Beschleunigung lings der Kugel, 


ein Maximum haben. 
Versuche mit Kugeln von 25 und 38 mm Durchmesser, aufgehingt in einem Rohr von 


146 mm lichter Weite, ergaben gute Ubereinstimmung mit der Theorie. Wellenbildung dieser 


Art kann die Belastung von Fillkérpersiiulen mit grossem Leervolumen beeinflussen. 


1. INTRODUCTION liquid leaving as a jet at the bottom. When the 


surrounding air is moved vertically upwards the 


Ix a previous paper [1] it was suggested that 


loading in packed towers is caused by the forma- velocity just above the liquid surface at A 


tion on the liquid film of standing waves supported 


by pressure gradients within the gas stream. Liquid saline F 


These pressure gradients are due to the accelera 


tion of the gas as it passes up through regions of aaa saa 

the packing where the flow area is contracting. ii] i | | 
The present paper deals with the formation Air \ if Flow 

of similar waves, but the pressure gradients were & < 

produced differently. Fig. 1 shows the arrange- a ® 

ment. A jet of liquid is put on to the top of a ji 


vertical string of spheres, and forms a symmetrical Fic. 1. Sphere covered by a liquid film in an upward 
film which runs down freely under 


gravity, the airstream, 
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will be less than at B. Consequently the pressure 
at A will be greater than at B, and the liquid 
film B will be 
tangential force caused by the pressure gradient. 
The to this 


pressure gradient is similar to that previously 


between 4 and subjected to a 


mechanism of wave formation duc 
described [1] ; 


fand B 


the foree of gravity must always be 


an element of the liquid film between 
in Fig. 1, 
balanced by the sum of the pressure gradient and 
viscous forces. The latter force is invers« ly pro 


portional to film thickness, and this will become 


very large at an airflow such that the gravity 
force is balanced entirely by the pressure vradient 
force. This balance produces the laree wave 
shown in Fig. 2 which has formed on the hemis 


pherical end of a 1-5 in. diam, cylinder covered by 


a liquid film; this arrangement was used to find 
the effect of altering the airtlow pattern above the 


lowest he musphe re. 


Expr riments were also made with vertical 
strings of spheres of diameters 1-0 and 1-5 in. 
The sph res and the evlinder wer SUS]x nded on 
the axis of a tube of 5-75 in. id., which is large 


enough to give a small wall effect. 

When the 
tube is small 1} the 
In this case the 
ball in the 


clearance between the balls and the 


wave formation ts different. 


gradient on the bottom 


pre ssuTe 


region 1B Fig. 1 is small compared 


with the pressure gradient in the region C. At ¢ 
the contracting area between the ball and the 
tube, as the air moves up, causes a pressure 
gradient, with wave formation just below the 


equator of each ball. This is thought to be the 
loading in towers of 
\ detailed treatment of this 
and will be published 
thought that 


tower, as 


mechanism of 


ked 
moderate voidage. 
case Is given elsewhere [2] 


With 


loading will begin at the bottom of the 


later. high voidages it is 


in the present experiments. 


2. EXPERIMENTAI 
Apparatus and procedure. 
A sectional scale drawing of the apparatus is 


shown in Fig. 3. 


Here A is the liquid distributor, 
B the string of balls, and C the tube of 5-75 in. i.d. 
The entry prece D was designed to have the same 
nozzle in B.S.1042 [3}. 


LS.A. The 


form as the 


liowkIns and 


if we consider the equilibrium of 
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i 
4 4 A 
4 
et A niet 
C 
4 
Liquic 
of packing sutiet inches 
ri rh apparatus (to seule) 
entry piece D fitted on top of a 1ft* packed 


tower, with a space of 1 ft between the bottom of 
D and the top of the packing. With these pre- 
cautions it was assumed that the air velocity 
bottom of the tube C. 


The liquid, after running as a film over the balls, 


would be uniform at the 


ran down a rod through the axis of the balls, and 
was collected by the small funnel and taken out 
by a small tube through the packing below. 

carried out with three 


Experiments were 


liquids ; water, water + 1 per cent Lissapol, and 
liquid paratlin. The water rate was measured by 
a rotameter, the liquid paraflin rate by timing a 
measured volume, and the air rate by a calibrated 


entry flume of the type described by Ower [4]. 


5. 
~ vO 
bes 7 

|| 


Wave forming on a film of liquid paraflin 


O-023 ml sec) due to the upward flow of au 


(velocity 17-5 ft sec) 


: 
‘ 
» 
ue 


q 
4 
: 
vo 
= 
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Countercurrent gas-liquid flow. 


With the liquid film running steadily, the air 


rate was gradually increased until a wave formed 


on the lowest hemisphere. 


3. THEORY 


The the outside of a sphere 


covered by a liquid film was analysed approxt- 


airflow round 


mately by using Stokes’ stream function, with 


various combinations of sources and sinks super- 


posed on a uniform stream, as described by 


Mitxr-Tuomson [5]. Typical air streamlines are 


shown in Fig. 4; here (a) shows the streamlines 


| / | 


(a) 


Inviscid motion (to scale). 


(b) 
(a) Flow 


Fic. 4. round a 


sphere, produced by a doublet in a uniform stream 


(b) Doublet and a point source in a uniform stream 


due to a doublet—that is, an equal source and 


sink at the same point combined with a uniform 


infinite stream. This combination of doublet and 


uniform stream gives the flow round a sphere. By 


appropriate choice of doublet strength and main 


stream velocity, one of the streamlines AA can 


be made to have nearly the same shape as the 


Over the lower half 


surface of the liquid film. 
of the sphere the actual air streamlines are similar 


to those shown in Fig. 4(a), and consequently the 


pressure gradients along the first part of the 


streamline AA will be similar to those which 


actually act on the liquid film. 
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On account of separation of the main air- 
stream. and the formation of a wake behind the 
sphere, an element of air starting in contact 
with the water surface flows along a path BB 
AA. the 


gradients over the lower half of the sphere is not 


rather than The change in pressure 
large and was estimated from the streamlines in 
Fig. ib). These streamlines are due to a comet 
dent source and doublet with a uniform stream. 

The stream function for the flow round a sphere 


(Fig. 4a) is: 
Ur? sin® @ 


—— | l 


a® 
Here U is the main stream velocity, r is the 
radius measured from the doublet, a ts the radius 
of the sphere, and @ is the angle measured from 
the axis. The streamlines of Fig. 4(b) are obtained 
by adding a term m cos @ to equation [1], m being 


the strength of the additional source at the origin. 


Boundary conditions 


(a) The streamlines 4A and BB are chosen to 
coincide with the surface of the film far below the 
sphere. This means that if R, is the radius of the 
rod plus the asymptotic film thickness, in Fig. 4 
the value of % along AA is UR,/2 and is 


(UR,?/2) + malong BB. 


(b) When @ = 90° the distance of the streamline 
from the axis must be R,, the radius of the sphere 


plus the liquid film thickness. For the streamline 


AA, the resulting equation ts = — a’, 
and for the streamline BB there is the same 
equation with 2m R/U subtracted from the 


right hand side. 


(c) For the streamlines in Fig. 4a the boundary 
(a) the 
completely _ but a third condition is required when 


conditions and (b) determine motion 
separation is considered as in Fig. 4b. Rather 
arbitrarily, it was decided to equate to R, the 
distance of the streamline BB from the axis at 
@ — 100°. The resulting equation is UR 

U(R,* — a® sin® 100°) — 2mR, (1 — cos 100°) and 
this, together with boundary condition (b), de- 
termines a and U/m. In this case a is nearly 


equal to the radius of the actual sphere. 


voL. } | | 
\ 
\ \ 
\ 
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Pressure gradient 


The radial and tangential velocities u and v 
along the streamlines AA and BB were 
calculated by substituting equation (1) into 
the expressions 0)/r sin @ and 
i (d¥/dr) rsin@. The pressure p at any 
point was then calculated from BrrNou u's 
theorem, and plotted as a function of @. The slope 
of the resulting graph is (dp /d@),. 


Condition for wave formation 

The airflow at which the pressure gradient forces 
entirely balance the gravity forces on the liquid 
film is calculated by putting : 


r\d@/, 

Here p, is the liquid density. This equation is 
derived by a balance on an element of the liquid 
film, and gives an upper limit to the airflow for 
stability. A similar upper limit calculated in the 
previous paper [1] for a vertical water film was 


about 25 per cent above the true criti al airtlow, 


Angie.) degree: 


Fic. 5. Functions for determining the critical airflow. 


Fig. 5 shows how equation (2) was satisfied, 
and gives sin @ and ()p/d0)/r p, gU* as functions 
of @. To satisfy (2), sin @ is multiplied by p,/p,U*, 
p, being the density of air. The value of U was 
adjusted to give contact between the curves in 
the region X Fig. 5. This is the critical value of U 
and the wave forms at the point of contact. In 
the case of the 1-0 in. sphere the angle is 29° and 
is indicated in Fig. 6 with the streamline BB 
which was assumed to coincide with the surface 
of the liquid film. In deriving (2), it has been 
assumed that the liquid film is very thin, that is, 
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that the streamline BB is nearly tangential to the 
ball. This is only approximately true at 29° in 
Fig. 6, so too much accuracy must not be ascribed 
to these calculations. 

The procedure when using the streamlines of 
Fig. 4a, not allowing for separation, is similar, 


and results from both methods are given in 
Table 1. 


\ 
29° 
x 
| fy 
\ 
| 
8 
Li 4 


Fic. 6. Air streamline round a 1 in. diam. ball. 


Liquid paraffin rate = 0-03 ml/sec. 


4. RESULTS 


Table 1 shows the results for two strings of 
spheres and for the cylinder with hemispherical 
ends. The latter was used to get an experimental 
measure of the effect of altering the airflow above 
the lowest hemisphere. Three liquids were used, 
water and liquid paraffin because of their widely 
differing viscosities, and water with 1 per cent of 
Lissapol because of its low surface tension, about 
30 dyn,cm, its other characteristics being the 
same as water. The liquid flow rates were 
adjusted to give roughly the same film thickness 
for each liquid. 

It was not possible to measure the angle at 
which the wave initially became large, because 
the wave broke into spray with water, and moved 
round the ball with liquid paraffin. Thus in Fig. 2, 
where the initial formation was at about 30° to 
the vertical, the wave as photographed is at 
about 45° to the vertical. 


B\ 
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Table 1. Wave formation on a sphere in an air stream 


Viscosity of the liquid paraflin=> 300 cP 


Theoretical 


Experimental 


Sphere Liquid - - — _— 

dia., Liquid rate, No separation Separation Col. 8 

in. mil /sec velocity, angle, velocity, argle lir velocity Col. 6 
ft/sec deg ft /sec deg ft /sec 


Liquid 0-03 233-3 
paraftin 
10 Water 35 
10 Water 1% 
Lissapol 
15 Liquid 
paraflin 0-03 


Water 


Water + 1% 
Lissapol 
Cylinder | 
dia. 
1-5 Liquid 
paraffin 0-023 
Water 


Water + 1% 
Liss 


32 24 29 17-9 O-74 
26 22-3 O-8T 

25°8 26 19-2 

24 17-3 0-61 


te 
- & 
— 


O-O4 


24 19-3 


5. Discussion 
Effect of separation behind the lowest sphere 


The results given in Table 1 from the two 
theoretical methods, the streamlines being those in 
Fig. 4(a) and 4(b), are in good agreement, showing 
that the effect of separation is small. Further 
evidence is provided by the excellent agreement 
between the results for the string of 1-5 in. 
spheres and the results for the 1-5 in. cylinder with 
hemispherical ends. In these two cases the airflow 
pattern above the lowest hemisphere must be 
different although the difference may not be great. 

These results are what would be expected, since 
the flow pattern round the nose of a streamlined 
object is not much influenced by the shape of its 
back. 

Another effect of separation is to ensure that the 
wave formation is only on the lowest sphere. 
Every other sphere is in the wake of the one 
below, and since a wake is a region of constant 
pressure, no wave is formed. 
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Interaction between the air and water streams 


In the previous paper [1] elaborate relaxation 


calculations were necessary to determine the 
critical airflow for wave formation, 0-115 ft*, sec. 
These calculations were necessary because the 


pressure gradients inducing the wave were in- 
fluenced by the shape of the wave and vice versa. 
When this interaction was ignored, as in the 
present paper, a value of 0-145 ft®/sec was ob- 
tained, the ratio of the two values being 0-79. 
This is of the same order as the ratio of the 
experimental to the theoretical velocities given 
in the last column of the Table. Therefore it 
seems likely that the difference between theory 
and experiment in the present case is largely due 


to interaction. 


Surface tension 

In the previous experiments [1], the critical 
airflow was raised about 5 per cent when Teepol 
was added to the water, and this was explained as 


| 
15 | 1-76 to = 24 22-6 to : 
8-7 | 23-8 
VOL. 30-2 
1957/58 
17°5 
22-6 to 
8-7 23-8 
1-5 
| | | | | | 
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being the result of eliminating travelling ripples. 


These ripples, in passing through the standing 


waves, were thought to cause premature in- 


stability. 
In the the 


eliminates travelling ripples, the critical velocity 


present case, although Lissapol 
is lowered about 13 per cent. A likely explanation 
is that the pressure gradients are quite sensitive 
to the shape of the liquid surface near Y, Fig. 6, 
and this shape is affects dl by surface te nmsion. This 
the 


paratlin 


explanation by 
the 


Lissapol, both liquids having a surface 


Is supported agreement 


between results for liquid and 
water 
tension around 30 dyn, cm. 

This sensitivity to the shape of the air stream- 
the Fig. 6, also makes the 


theoretical results unreliable. they 


lines in region Y, 


since are not 


based on the shapes of the actual liquid surfaces. 
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NOTATION 


air stream velocity 
radius of rod plus film 
radius of sphere plus film 
radius of sphere 
strength of source 
pressure on the liquid film 
radial distance 
radial velocity 
tangential velocity 
angle between the radius vector and the axis 
density 
stream function for the airflow. 
S ripts 
l liquid 


a air 
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Vapour-liquid equilibrium still for low pressures 


LYDERSEN 


The Technical University of Norway, Chemical Engineering Laboratory, Trondheim, Ni rway 


August 1957) 


and Eci. Hammer 


\ new still for vapour-—liquid equilibrium measurements at pressures down to 1 mm Hg 


is described. The whole still is agitated in order to give vaporization from a frequently renewed 


falling film, and also to give good mixing of the boiler liquid and the condensate. 


Possible sources of error are reviewed and tests were made to check on the different types 
of error. The enrichment obtained in the still appears to be approximately 98 per cent of the 


enrichment by complete equilibrium. 


Vapour-liquid equilibrium measurements 
di-n-butylsebacate at 1 mm Hg are given, showing practically ideal behaviour of the mixture. 


Résumé 


with a mixture of di-n-butylphthalate and 


Les auteurs décrivent un nouvel ébulliométre pour mesurer des pressions inférieures 


a ilmm Hg. Tout lappareil est agité de fagon & provoquer l’ébullition 4 partir d'un film liquide 


constamment renouvelé, et aussi pour favoriser un bon mélange du liquide bouillant avec le 


condensat. 


Les auteurs passent en revue les diverses sources d’erreurs et font des essais pour discerner 
ces diverses causes d’erreur. L'enrichissement obtenu dans l'ébulliométre semble voisiner 98 pour 


cent de lenrichissement obtenu a l’équilibre. 


Les mesures de l’équilibre liquide vapeur, faits avec un mélange de di-n-butylphtalate et 
de di-n-butylsebacate 4 1 mm Hg, montrent le comportement pratiquement idéal du mélange. 


Zusammenfassung— Es wird eine neue Apparatur zur Messung des Dampf-Fliissigkeits—Gleich- 
gewichtes bei Driicken bis herunter zu 1 Torr beschrieben. Die gesamte Apparatur ist geriihrt, 
um eine Verdampfung aus einem hiiufig erneuerten Rieselfilm zu erzielen und um auch eine gute 
Durchmischung der siedenden Fliissigkeit und des Kondensats zu erreichen. 

Die mdglichen Fehlerquellen werden diskutiert und an Probeliufen die verschiedenen 
Fehlerarten untersucht, Die mit der Apparatur erzielte Anreicherung scheint ungefihr 98 °% 
von jener zu betragen, die bei vollstindigem Gleichgewicht zu erreichen wiire. 

Messungen des Dampf—Flissigkeits-Gleichgewichts einer Mischung von Di-n-Butylphthalat 
und Di-n-Butylsebacat bei 1 Torr werden mitgeteilt und zeigen das praktisch ideale Verhalten 


des Gemisches. 


INTRODUCTION 
THe growing use of vacuum distillation has 
increased the importance of knowing low-pressure 
vapour-liquid equilibria. The experimental deter- 
mination of equilibria at low pressures involve 
a number of difliculties not encountered at 
elevated pressures, primarily because only minor 
pressure drops can be tolerated, but also to some 
extent due to the tendency of bumping. 

Fow er [1] gives a description of the vapour 
liquid equilibrium stills published in the literature 
up to 1948 and RrpGway (9) gives a correspond- 
ing review for the years 1948-55. Recent modi- 
fications of the stills are described by several 
authors ([8], [4], [5] and [10)}). 


241 


The majority of the stills do not seem adaptable 
to pressures below a few mm of Hg, as the bubble 
formation and growth take place under an 
appreciable liquid head. At a pressure of 1 mm Hg 
in the vapour space, the pressure 2 cm under the 
liquid surface will be 2 to 3mm Hg, depending 
on the density of the liquid. 

The aim of the present authors was to build a 
still suitable for equilibrium measurements at 
pressures down to ImmUHg absolute. For 
dynamic stills only two main principles seem 
worth considering : 

The use of an inert carrier gas, whereby a 
total pressure is maintained that is high compared 
to the partial pressure of the compounds in the 
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liquid mixture, thus making it possible to use a 
conventional type of still. 

Vaporization from a liquid surface that is 
renewed fast enough to avoid any appreciable 
dilution with respect to the more volatile 
compound, 

In the present work, the design was based on 
the latter principle, mainly because, firstly, the 
vapour pressure of the mixture can be measured 
directly and, secondly, because the total condensa- 
tion of pure vapours is easier to obtain than the 
condensation from an inert, where fractional 
condensation will cause error. 


L. Lypersen and Eom HamMer 


DESCRIPTION OF THE STILL 

The still to be described in this paper is shown 
in Figs. 1 and 2. The vapours formed in the 
boiler A, Fig. 2, are condensed in the parts D 
and E, and the condensate flows back to the 
boiler through the capillary G. Heat flow 
through the inner half of the insulation of the 
off-take pipe C is prevented by electric heating 
of the aluminium cylinders H, and H,. The heat 
input is manually adjusted according to the 
readings of the thermocouples 1-7. 

The still is connected to a McLeod gauge and a 
rotary vacuum pump protected by the dry ice 
traps K, and K, as shown schematically on the 
right-hand side of Fig. 2. An air-bleed valve N 
and a Cartesian manostate V makes it possible 
to adjust the vacuum. 

The top of the still is fastened to a universal 
joint U. During operation, the still is agitated 
as the shaft Y is rotated in a circle with radius r. 
This gives agitation of both the boiler liquid and 
the condensate in the sphere F. The capillaries 
for sampling and liquid return serve as baffles. 

The necessary operations are adjustment of 
the vacuum and the heating elements in the 
insulation, and also checking temperatures and 
pressure to see that stationary conditions are 
maintained. At the end of the run, the vacuum 
is released. This instantaneously stops the boiling 
and forces boiler liquid and condensate into the 
sample flasks R, and R,. 

The design shown on the Fig. 1 and 2 has some 
resemblance to the still used by Witi1aMs [11], 
which again was a modification of an OrHMER 
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Fig. 2. A drawing of the still; the connections at the 
lower right are shown schematically, Legend : A = 6 em 
dia. boiler with heating mantle B, C = vapour off-take 
pipe leading to aircooled condenser D and the water- 
cooled condenser E, F = condensate receiver with over- 
flow capillary G; H, and Hg, = electrically heated alu- 
minium cylinders ; K, and Ky = dry ice traps protecting 
the McLeod gauge M and the vacuum pump P; N is an 
air-bleed valve, T a surge tank and V a Cartesian 
manostate ; R, and R, = sample flasks for condensate 
and boiler liquid, and Y = shaft, agitated with radius 
r by means of the rotating disk Z. 


[7] still. Wititams, however, used a filling of 
helices in the boiling chamber to avoid bumping, 
while Lancpon and Keyes [6], Perry and 
[8] and later Sanpers and Acciarri 
[10] made use of mechanical stirrers. 

The type of agitation used in the present still 
causes the bulk of the liquid to rotate along the 
side of the boiler, leaving a falling liquid film 
on the rest of the boiler surface. This ensures 
vaporization from a liquid film that covers the 
major part of the boiler surface which at any 
moment is not covered by the bulk liquid. This 
additional surface is considerably larger than the 
surface over the bulk liquid. 
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Vapour-liquid equilibrium still for low pressures 


ESTIMATION OF ERRORS 


Inadequate renewal of the vaporizing surface 
might cause the vaporizing film to be leaner in 
the more volatile compound. This type of error 
will decrease both with decreased vaporization 
rate and with increased speed of the agitating 
disk. Fig. 3 shows the enrichment in the volatile 
compound obtained from a 50 mole per cent boiler 
mixture of di-n-butylphthalate and_ di-n- 
butylsebacate as a function of the boiling rate at 
three different speeds of the agitation. With a 
radius r (Fig. 8) of 9mm a maximum of 
260 rev/min and a minimum distillation rate of 
approximately 0-35 kg/hr was chosen for easy 
operation of the still. This boiling rate corres- 
ponds to vaporization of the boiler content 10 
times an hour and renewal of the condensate 
sample 150 to 200 times an hour. 


Fic. 8. The enrichment of di-n-butylphthalate ; y-z, 
obtained from a 50 mole per cent boiler mixture of 
di-n-butylphthalate and di-n-butylsebacate at 1 mm Hg 
as a function of distillation rate W kg/hr by three different 
speeds of rotation n of the disk Z shown in Fig. 2. 


Extrapolation in Fig. 3 to zero distillation rate 
indicates that the enrichment obtained at the 
chosen approximately 
98 per cent of the enrichment at equilibrium 
conditions. 

Rectification effect can occur by improper 
heating of the insulation of the boiler off-take 


conditions should be 


pipe. However, changing the temperature drop 
outwards through the inner half of the insulation 
from 0 to 7°C gave no measurable change in the 
composition of the condensate. 

Total vaporization of boiler liquid from the 
boiler off-take walls may result from splashing 
or entrainment from the boiler in connection 


with superheating of the pipe C shown in Fig. 2. 
Tests with a temperature drop inwards through 
the insulation of 7°C reduced the enrichment in 
di-n-butylphthalate to only 82 per cent of the 
enrichment without superheating. Superheating 
only the upper half of the pipe C, however, 
gave no measurable change in condensate com- 
position. This indicates that droplets of boiler 
liquid only reach the lower part of the boiler 
off-take pipe, and that safe operation of the still 
requires a slight temperature drop outwards 
through the lower half of the insulation. 

Liquid entrainment carried over to the condenser 
does not seem probable as long as no vaporization 
of entrainment could be found from the upper 
part of the off-take pipe. Moreover, the correction 
derived from Fig. 3 should include the influence 
of any liquid carry-over. 

The condensate sample may not be representative. 
Condensate draining from the condenser after the 
boiling stops may have a composition deviating 
from the average composition of the condensate. 
The equilibrium measurements reported in the 
next section indicates no systematic error. Hence, 
any serious error of this type does not seem 
reasonable. 

Flashing of the low boiling condensate could 
occur as the condensate returns to the boiler. 
Substantial subcooling of the condensate and 
returning it under the liquid surface is believed 
to make errors of this type negligible. 

Vaporization under a liquid head seems not to 
occur; visible bubble formation disappears 
immediately when the agitation is started. 

Contamination with stop-cock grease is unlikely 
to occur, as no liquid comes in direct contact 
with the grease used on the conical joints of the 
sample flasks. Teflon gaskets were used for the 
vacuum and gauge connections. 

Test runs were carried out to check the liquid 
filling necessary for correct operation. Only 
8 cm® of liquid in the boiler gave an appreciable 
error, probably due to the fact that the bulk 
of the boiler liquid did not wash over the upper 
parts of the boiler walls. Boiler fillings from 
82 to 55cm*, however, gave consistent results. 
In addition there is a condensate hold-up of 
2-8 cm’. 
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VarpourR-Liquip EQuILIBRIA 
MEASUREMENTS 
The new still was used to determine the vapour- 
liquid equilibria for mixtures of di-n-butyphthalate 
and di-n-butylsebacate at a pressure of 1 mm Hg. 
This mixture is satisfactory for the calibration 
of low-pressure distillation columns with less 
The 


easily ay ailable. and the mixtures can be analysed 


than 6 theoretical stages. chemicals are 


by measuring the refractive index. 
Also, it to that 
mixture should show an almost ideal behaviour. 


is reasonable believe this 
Previous vapour pressure measurements [2] gave 
a relative volatility « as given by the equation : 


26-6 84800 


— —- +- 
7 T 


log « (1) 
At I mm Hg, the boiling point of the two com- 
pounds in question is 148-8 and 164-1°C respec- 
tively, and therefore « calculated from equation 
2-31 and 2-59. 


The materials used were the same as for the 


(1) varies between 


vapour pressure measurements [2]. The results 


L. Lypersen and Kou. Hammer 


of the equilibrium measurements are reported 


in Table 1 and plotted in Fig. 4. The mole 


1-0; 
y 


Fic. 4. 
mixture of di-n-butylphthalate and di-n-butylsebacate 


The vapour-liquid equilibrium diagram for a 
at 1mm Hg. Here curve a gives the mole fraction y of 
the phthalate in the vapour as a function of the mole 
fraction 2 in the liquid as calculated for an ideal mixture. 
The open circles represent the experimental data. Curve 
b gives the temperature as calculated for an ideal mixture, 
the black dots representing the thermocouple readings. 


Table 1. Vapour-liquid equilibria for di-n-butylphthalate di-n-butylsebacate miatures at 1mm Hg. 


3 4 


Vapour temperature (°C) 


Mole fraction of di n-butylphalate 


Pressure 


Run No, in vapour 


(mm Hg) in liquid 


= lideal 


Yexp Videal 


0-139 
0-269 
0-272 
0-386 
0-474 
0-495 
0-579 
0-710 
O-718 
0-792 
0-850 
0-874 
0-880 
0-940 
0-951 
0-954 


162-6 
161-5 
161-4 
159-4 
158-1 
158-4 
156-6 
1551 
154-3 
153-2 
152-4 
151°8 
151-1 
151-0 
150-0 
150-4 


162 
161: 
161-2 
159-7 
158-6 
158-3 
157-1 
154-9 
154-8 
153-4 
152° 
151-7 
150-3 
150-0 
150-0 


0-143 
0.281 

0-279 
0-388 
0-477 
0-496 
0-582 
0-712 
0-716 
0-788 
0-842 
0-867 
0-874 
0-038 
0-947 
0-950 


0-065 
0-136 
0-138 
O211 

0-276 
0-293 
0-366 
0-505 
0-514 
0-611 

0-608 
0-740 
0-750 
0-864 
0-888 
0-804 


1-000 
1-005 
1-000 
1-000 
0-995 
1407 
0-998 
1-005 
1-000 
0-998 
1-000 
1-000 
1-000 
1-001 
1-000 
1-000 
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fraction of phthalate in the vapour, y,,, was 
calculated using the equation : 


Vexp — = (y — 2)/0-98 (2) 


where y is the mole fraction actually measured 
in the condensate, and the factor 0-98 corrects 
for the deviation from equilibrium. 

The mole fraction in Table is 
calculated from the relative volatility given by 
equation (1) assuming ideal behaviour. The 
temperature f,.., in Table 1 is the boiling 
point calculated for an ideal mixture under a 
total pressure of 1 mm Hg. 

Experimental vapour compositions and ex- 
perimental boiling points are plotted in Fig. 4 
fraction of di-n- 
The curves 


as a function of the mole 
butylphthalate in the boiler liquid, 2. 
a and b in Fig. 4 correspond to values calculated 
for ideal mixture behaviour. The experimental 
points correspond to the curves for ideal behaviour 
within the experimental error. Hence, vapour 
liquid equilibria data for the mixture in question 
can readily be obtained from equation (1) com- 
bined with the previously published [2] vapour- 


pressure data. 


Vapour-liquid equilibrium still for low pressures 
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NOTATION 

a = relative volatility of di-n-butylphthalate and 
di-n-butylsebacate 

n = number of rev/min for the disk Z shown in 
Fig. 2 

r = mole fraction of di-n-butylphthalate in the 
boiler liquid 

y = mole fraction of di-n-butylphthalate found in 

the condensate 

-mole fraction of phthalate in the vapour 


Yexp 
calculated from equation (2) 
Yidea) = Mole fraction of phthalate in the vapour as 
calculated for an ideal mixture 
texp = measured temperature of the saturated vapour, 
lideal = temperature as calculated for an ideal mixture, 


T = absolute temperature, °K 
W = rate of distillation, kg/hr 
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Vapour-liquid equilibria in the ternary system 
ethyl acetate—benzene—cyclohexane 


K. C. Cuao and O. A. Houcen 


Chemical Engineering Department, University of Wisconsin 
( Received 1 August 1957) 


Abstract—The vapour-liquid equilibrium compositions and temperatures covering the entire 
composition range were determined experimentally under isobaric conditions of 760 mm Hg 
pressure for the ternary system ethyl acetate — benzene — cyclohexane and its three constituent 
binaries. 

A continuous recirculating equilibrium still combining the best features of Jones and Gillespie 
stills was used. 

Values of activity coefficients and their ratios were correlated by Redlich—Kister equations. 
\ modification of these equations was introduced to satisfy the Gibbs—Duhem equation at 
constant pressure. 


Résumé—Les auteurs ont étudié expérimentalement, sous 760 mm Hg, les courbes de rosée et 
d'ébullition du systéme ternaire acétate d'éthyle—benzéne-cyclohexane et des trois mélanges 


binaires. 

Iis ont utilisé un ébulliométre dynamique présentant les principales caractéristiques des 
ébulliométres de Jones et Gillespie. 

Ils ont déterminé les valeurs des coefficients d'activité et leurs rapports par les équations de 
Redlich-Kister. 

Ils ont introduit une modification de ces équations pour satisfaire A la relation de Gibbs- 
Duhem. 


Zusammenfassung—Zusammensetzungen und Temperaturen des Dampf-fliissig—Gleich- 
gewichts dber den ganzen Bereich der Zusammensetzungen wurden experimentell unter isobaren 
Bedingungen bei 760 Torr fiir das ternire System Athylazetat-Benzol-Cyclohexan und die drei 
darin enthaltenen Zweistoffsysteme bestimmt. 

Verwendet wurde eine kontinuierliche Umlauf-Gleichgewichtsapparatur unter Ausnutzung 
der besten Vorschlige von Jones und Gillespie. 

Werte der Aktivitatskoeffizienten und ihrer Quotienten wurden durch Gleichungen nach 
Redlich und Kister wiedergegeben. Diese Gleichungen waren abgeiindert, um die Gibbs-Duhem— 
Gleichung bei konstantem Druck zu befriedigen. 

OBJECTIVES behaviours of multicomponent systems and those 


Most of the published data on vapour-liquid of lower-order constituent systems. 


equilibria pertain to binary systems. Data This investigation has the following objectives: 


covering the entire composition range on relatively 
few systems have been published. From an 
industrial standpoint experimental data on 
ternary systems are essential because of their 
frequent incidence in chemical processing and 
their growing importance in azeotropic distillation 
and solvent extraction. From a_ theoretical 
standpoint such investigations are also important 
in establishing useful relations between the 


To determine the vapour—liquid composition— 
temperature relations at 760 mm Hg pressure of 
the ternary system ethyl acetate—benzene- 
cyclohexane and its three binaries; to correlate 
the experimental data of activity coefficients of 
the ternary system and its three binaries in terms 
of composition; and to compare the correlation 
constants in equations for isobaric vaporization 
in conventional and corrected forms. 
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The particular ternary system ethyl acetate— 
benzene—cyclohexane was selected because of its 
possible industrial importance in the separation 
of the benzene—cyclohexane azeotrope by dis- 
tillation with ethyl acetate. The mixture of 
benzene and cyclohexane is an industrial product 
resulting from the catalytic hydrogenation of 
benzene to form cyclohexane. 

The addition of ethyl acetate to the binary 
system benzene—cyclohexane in fractional dis- 
tillation results in the separation of benzene from 
the binary azeotrope of cyclohexane-ethyl 
acetate. The separation of ethyl acetate from this 
latter azeotrope may then be made by solvent 
extraction with water. 


Previous LITERATURE 

The three binary systems reported in_ this 
research have been investigated by Tao [19] at 
pressures of 1, 10, 20 and 30 atm. The benzene— 
cyclohexane system has been investigated by 
Scarcuarp, Woop and Mocne. [16] at 39-99° 
and 69-98°C and by Ricuarps and HARGREAVES 
[14] at 759mm Hg. Furfuraldehyde and methyl 
cellosolve in ternary systems with benzene 
cyclohexane have been studied by THorNnTON and 
Garner [20], [21]. The ternary system selected 
for this investigation has not been previously 


reported in the literature. 


EQUILIBRIUM STILL 


The continuous equilibrium still used in this 
investigation (Figs. 1 and 2) is essentially a 
combination of the Jones [9] and GiLLespre [5] 
types. 

With reference to Fig. 1, contact between 
vapour and liquid is secured in the mixing 
chamber A from which vapour with entrained 
liquid rises to impinge on the thermowell W, 
located in the separator D. The liquid falls to 
the trap E, of 23 ml capacity. Liquid overflows 
from this trap to the mixing chamber thus 
completing the circuit for the liquid. Vapour 
escaping from the separator is totally condensed 
in condenser F and collected in the condensate 
trap G from which the condensate flows at a 
closely controlled rate into the vaporizer K. 
Upon complete vaporization the vapour bubbles 
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Fic. 1. Equilibrium still. 


through the liquid in the mixing chamber 4, 
completing its circuit. 

To obtain significant data it is essential to 
circulate the liquid under adiabatic conditions. 
For this purpose an insulating jacket encloses 
the liquid trap and the mixing chamber. 
The innermost part of the jacket 
of a_ sheet-metal box L with dimensions 
2} x 4} x 10in. Two layers of asbestos fabric M 
are wrapped around this box with a chromel 
ribbon heater H, placed in between the layers. 
An outer box made of plywood jin. thick 
encloses the assembly (34 x 6 x 10}in.). The 
space inside the sheet metal box is filled with 


consists 


Santocel insulation. 

The vapour-liquid tube C that extends above 
the jacket is wrapped with glass wool and wound 
with a heating wire H, to permit adiabatic per- 
formance. 

Auxiliary iron—constantan thermocouples for 
purposes of adiabatic control are attached to the 
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wall of the mixing chamber, of the residue liquid 
trap and to the different sides of the sheet-metal 
box. 

The equipment is heated electrically with four 
chromel ribbon and wire heaters controlled by 
three Variacs. Heater 1 consists of chromel-A 
ribbon (4 x 0-04in., resistance 2-44 Q ft, the 
total resistance being 160 2) wound uniformly on 
all four sides of the adiabatic jacket which encloses 
the entire liquid recirculating system. Heater 2 
A wire (size 30: 6-71 Q ft. 
the total resistance being 68 2) supplies heat 


consisting of chromel- 
to the vaporizer. Heater 3 supplic s heat to the 
bottom of the adiabatic jacket. This heater 
consists of chromel-A wire having a total re- 
sistance of 20 2 as in heater 2; the wire is 
enclosed between two layers of asbestos fabric as 
in heater 1. Heater 4 maintains adiabatic condi- 
tions in the vapour-liquid tube C which extends 
above the box enclosure. 

To promote even heating a sheet of aluminium 
foil surrounds the vaporizer; _ the foil is separated 
from the heating coil by a layer of glass wool. 
A closely wound coil of 22 gauge copper wire fits 
tightly against the inner wall of the vaporizer to 
act as a liquid distributor. This coil induces a 
liquid film to form and wet the entire inner 
surface of the tube. 

To maintain the still at constant pressure, a 
regulating device is attached to the top opening of 
the condensate trap of the equilibrium. still 
through an auxiliary condenser. This control 
is effected by bleeding compre ssed air through a 
nozzle immersed in water, the depth of immer- 
sion being adjusted to compensate for variations 
With this device 


the pressure in the still may be kept constant 


in the atmospheric pressure. 


within an accuracy of 1-5 mm of water. 

The equilibrium temperature was measured by 
a calibrated copper constantan thermocouple T', 
located in thermowell W, in the separator; the 
thermowell being filled with liquid silicone for 
thermal contact. The 
readings was within 0-1°C. Readings of thermo- 


improved accuracy of 
couple 7, located in the mixing chamber A were 
made for comparison with those of couple T',. 

All auxiliary thermocouples were iron con- 
stantan recording to an accuracy of 1 F. 
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OPERATION OF THE EQUILIBRIUM STIL! 


In starting operation the adiabatic jacket and 
the tube to the separator are heated to the antici- 
pated temperature, about 2 hr being required for 
preheating. Cooling water is admitted to the 
total condenser, and heat applied to the vaporizer. 
A voltage of 35 to 40 V is applied to heater 2 for 
about 5 min before admitting feed to the still. 
The auxiliary condenser is removed from the still 
and the stopeock leading to the vaporizer is 
slightly opened. Then a charge of known composi- 
tion is fed to the condensate trap through a 
connection leading to the auxiliary condenser. 
A mixture of liquid and vapour passes into the 
about 68 ml of total liquid 
charge is required. The liquid should be supplied 


mixing chamber ; 


gradually so that boiling may be continuous in 
the vaporizer; this insures saturation of the 
liquid in the liquid trap as soon as the charging 
operation is completed. Feeding usually takes 
about 15 min, after which the stopcock S, leading 
to the vaporizer A is partly closed to its normal 
operating position. The auxiliary condenser is 
restored and the cooling water circulated. <A 
constant pressure of 760mm Hg is applied to 


the still. 


During operation adiabatic conditions are 
maintained by adjusting the adiabatic heaters 
H,, Hs, Hy so as to equalize the temperature 
readings of thermocouples 7’, and 7, and those on 
the adiabatic jacket. Thermocouples 7, and T, 
show no disagreement after a short period of 
operaton. Temperatures on the different parts of 
the adiabatic jacket may vary by 2°F, a range 
which is adjusted to include the temperatures of 
thermocouples 7, and 7; thus the temperature 
difference across the Santocel insulation is in 
every case less than 2°F. 

Heat to the vaporizer is controlled so as to 
cause the last trace of condensate to disappear 
just at the end of the vaporizer, a condition 
which can be observed visually. The sensitivity 
of this observation is enhanced by the V-shape of 
the vaporizer permitting traces of liquid to be 
easily seen in the bottom of the V. 

Equilibrium conditions were usually reached 
within 1 hr of operation. It was customary to 
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hic, 2. Photograph of equilibrium still. 


Fic. 3. Precision pycnometer. 
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Vapour- liquid equilibria in the ternary system ethyl acetate—benzene—cyclohexane 


operate 90min before taking equilibrium 
readings. 

At the end of each run, the stopcock S, leading 
to the vaporizer is closed and the vaporizer 


off. 
condensate trap and the liquid trap are opened 


heater turned The stopeocks below the 
to sampling bottles after purging the line with 
8ml of the liquid. 


withdrawn the entire content of the still is emptied 


When samples have been 


in preparation for the next run. 

No difficulty was encountered in securing leak- 
proof operation of the two stopeocks 8, and S, of 
the condensate trap; these stopcocks were not 
greased. However, difficulty was experienced 
with the two stopcocks S, and S, below the mixing 
found 


chamber and the liquid trap. It was 


necessary to use high vacuum grade cocks to 
obtain good service ; these required lubrication. 
A thin layer of silicone lubricant was applied, 
but it was necessary to assure that no contamina- 
tion resulted from this source. To verify the 
absence of contamination the three pure liquids 
used in this study were tested separately. The 
first 5 ml of liquid withdrawn from the sampling 
port were collected, but only in the case of 
benzene was the first 5 ml sample found different 
from the liquid in the trap. Further samples 
showed exactly the same refractive index as the 
original feed. Contamination of the first 5 ml in 
the case of benzene corresponds to a change of 
0-0001 in 


connection to the greased stopcocks prevented 


refractive index. The capillary tube 


contamination by diffusion. To insure no con- 
tamination 8 ml of liquid were purged from the 
line before liquid samples were taken for analysis. 


ESTABLISHMENT OF EQUILIBRIUM 


An inherent advantage of the vapour recir- 
culation still is that the two phases retain their 
separate compositions and flow rates, with no 
net vaporization or condensation. Thus a 
steady state attained under conditions of recir- 
culation is an unmistakable sign of equilibrium. 
The following five conditions of operation were 
essential to assure the establishment of equili- 
brium conditions. 

(1) One manifestation of steady state is the con- 


stancy of the boiling point of the liquid as 


measured by thermocouple 7',. Variation in this 
reading usually disappears after about 30 min 
operation. 

(2) Constancy of composition of the vapour and 
liquid streams throughout is another indication of 
steady state and equilibrium. Since the amount 
the still 
maintained within close limits it is not possible 


of material circulating in must be 
to withdraw samples sufficiently large for analysis 
and to continue operating at the same time. The 
volume of liquid in the still is critical for proper 
operation of the vapour lift, and the amount of 
condensate is critical in furnishing the proper 
hydrostatic head to force the liquid into the 
vaporizer. A series of runs were made on about 
the same mixture of benzene and ethyl acetate in 
order to determine variation in composition of 
liquid in the still with residence time. 

Mixtures of benzene and ethyl acetate con- 
taining about 30 mole per cent of ethyl acetate 
were charged to the still. Successive runs of 
different 


results are given in Table 1. 


Table 1. 


time intervals were made and the 


Length of operation “4 
(min) 


v4 | Run number 


AB-19 
AB-17 
AB-16 
AB-18 


30 0.308 
0-303 
0-274 
0-288 


0-025 
0-028 
0-029 
0-029 


Even though there is no trend in 2 , values, the 
trend in y, — 24 is evident. Steady state condi- 
tions were obtained in from 40 to 60 min. It was 
concluded that 90 min of operation gave an ample 
margin of safety. 

(3) That equilibrium was attained at the end of 
the runs was further confirmed by withdrawing 
from the chamber. In all 


acceptable runs the samples showed the same 


samples mixing 
composition as those taken from the liquid trap. 
Identity of composition in the liquid trap and in 
the mixing chamber indicates that equilibrium 
was established. 

(4) In all instances constancy of the liquid level 
in the condensate trap gives an unmistakable 
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sign of thermal balance of the liquid recirculation 
system. This is one of the prerequisites to equili- 
brium and must be carefully watched throughout 
operation. In all acceptable runs change in the 
amount of condensate in the trap was scarcely 
noticeable, and never exceeded 0-2ml. Even 
under the most unfavourable conditions where 
this change would take place instantly the 
maximum error that could be made with the-most 
sensitive of the three systems would be 0-1 percent. 

(5). The vapour lift action in this still necessarily 
induces a temperature change due to variation in 
hydrostatic pressure with elevation. This pressure 
reduction is accompanied by vaporization of a 
portion of the liquid, and this would continuously 
increase the inventory in the condensate trap and 
deplete the inventory in the liquid trap. The 
change in saturation temperature due to hydro- 
static head in the still is slight. Thus temperature 
measurements at the bottom and top of the 
vapour lift, using thermocouples T, and Tp, 
showed a difference of 0-1 to 0-2°C between the 
two locations under steady conditions. That 
this hydrostatic effect is not appreciable, but 
comparable to that of imperfect adiabatic 
conditions around the liquid recirculation 
system, Is conclusively shown by the constancy of 
the amount of material in the two separate 
streams. 

The equilibrium pressure was maintained at a 
standard atmospheric pressure for all runs. Any 
discrepancy between the prevailing atmospheric 
pressure and the standard pressure was compen- 
sated for by applying compressed air to the 
system. That the pressure applied to the system 
though the auxiliary condenser actually represents 
the equilibrium pressure of the mixture depends 
on the negligible pressure difference for the 
vapour flow in the separator and in the condenser. 
This pressure drop was calculated not to exceed 
0-0081 cm of water. 


ANALYTICAL PROCEDURE 
Analyses of the three binary systems are 
readily carried out by measurements of refractive 
index and specific gravity. Measurements of these 
physical properties against mixtures of known 
composition were made. 
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For measuring refractive indices a Baush and 
Lomb Abbé refractometer was used. All measure- 
ments were performed at 25 + 0-1°C, the preci- 
sion of refractive index readings being + 0-000). 

For specific gravity measurements a modified 
Sprengel type pyenometer was used (Fig. 3). 
This weighed 17 g and held 10 ml of liquid in a 
single bulb with capillary connection tubing of 
0-7 mm bore. The reproducibility of the specific 
gravity measurement was +. 0-0001 at 25°C. 
The compositions of binary systems determined 
by these two measurements did not differ by 
more than 0-1 mole per cent. 

The physical properties of the standard binary 
mixtures are shown in Tables 1, 2, and 3, and in 
Figs. 4 and 5. 
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x 
VA 

N | 

\ | 

0-80 | 


O-2 O-4 06 O08 
Mole fraction of component 1 


Fic. 4. Specific gravities of binary mixtures. 


For calibration purposes specific gravities and 
refractive indices of ternary mixtures of known 
composition were measured over the entire com- 
position range at intervals of 10 mole per cent for 
a total of 66 mixtures; the calibration data are 
shown in Table 4. By linear interpolation 
isometric lines were constructed representing loci 
of constant values for each physical property. 
(Fig. 6). The isometric lines all merge into ter- 
minals on the edge of the triangular diagram 
coincident with data for the binary mixtures. 
Linear interpolation between the contour lines 
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Cyclohexane (1)-Ethy! acetate 
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Index of refraction 
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Fic. 5. Refractive indices of binary mixtures. 
gives the composition of an unknown sample 
once its refractive index and specific gravity 
were determined. 


0-87 0-89 
Ethyl! acetate 


0-8) 0-83 0°85 
Specific gravities 


0-79 
Cycionhexane 


Fic. 6. Refractive indices and specific gravities of the 
ternary mixtures. 


The precision of determining composition is 
+ 0-1 mole per cent for the binary mixtures and 
0-2 mole per cent for the ternary mixture. 


PREPARATION OF MATERIALS 


The ethyl acetate of 99-9 per cent purity was 
The chief 


impurity, water, was removed by storage over 


obtained from Merck and Co., Ine. 
silica gel. 

Benzene, of 99-9 per cent purity, was also 
obtained from Merck and Co. Further purifica- 
tion was obtained by fractional crystallization 
repeated three times. 

Cyclohexane of Research Grade (Lot Number 
435), of guaranteed purity 99-94 per cent was 
obtained from the Phillips Petroleum Co., and 
used directly. 

The physical properties of the purified reagents 
as measured are tabulated in Table 2 for com- 
parison with literature values. 


Table 2. 


Ethyl 
acetate 


( yelo 
Benzene hexane 
Refractive 
index at 25 
1-4235 
1-42354 [3] 


1-3698 
1-37005 [7] 


1-4980 
149790 [2] 


Experimental 
Literature 
Specific 
25) 
Experimental 
Density (25 ) 


Experimental 


gravity (25 
0-8970 O-8762 0-7759 
0-8736 0-736 

0-87362 [10] 0-7739 [4] 


O-8944 
Literature 0-89446 [8] 
Boiling point “¢ 

80-7 


80-7 [17] 


Experimental 


2 80-1 


Literature 7 [17] 80-1 [17] 


AZEOTROPES 


Because of the very close compositions of liquid 
and vapour in the equilibrium mixtures of ethy! 
acetate and benzene at high concentrations of 
ethyl acetate there is some uncertainty about the 
existence of azeotropism for this system. RyLAND 
[15] in 1899 reported the absence of azeotropism. 
Lecat [11] in 1918 reported an azeotropic com- 
position of 93 mole per cent of ethyl acetate, but 
the reliability of this finding was questioned by 
the same author. SwieTostawskt [18] in 1932 
confirmed Ry.anp’s finding. Tao [19] in 1952 
reported an azeotropic composition of 91 mole 
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per cent of ethyl acetate. To resolve the un- 
certainty, both Tao's data and the data of this 
study were plotted in terms of y, y, against 2 4. 
On the same figure was plotted the ratio P, P, 
against a,, the experimental temperature being 
used for calculating the vapour pressure ratio. 
Azeotropism would be indicated by intersection 
of the two lines. The two sets of data agree up to 
= 0-92. The temperature measurements also 
agree to within 0-1°C in both 
Judging by the trend of the data it is unlikely that 


investigations. 


an azeotrope exists. The one-constant equation 
giving a maximum deviation of 0-1 mole per cent 
for all the experimental measurements supports 
this observation in not giving any azeotropic 
composition. 

For the binary 
the reported values of azeotropic composition 


system benzene cyclohexane 
and temperature at 1 atm pressure are shown in 
Table 3. 

Table 3 


Benzene 


(mole Temperature 


per cent) Cc 


Lecat [11] 
RicHarps and 
HARGREAVES 
MARSCHNER and 
Cropper [12] 
THORNTON and GARNER 
Tao [19) 
This investigation 


For the binary ethyl 
acetate the reported values of the azeotropic 


system cyclohexane 
composition and temperature at 1 atm pressure 
are given in Table 4. 
Table 4. 

Acetate 


(mole 
per cent) 


Temperature 


Lecat [11) 53 
Tao [19) 53-3 
This investigation 54-9 


No azeotrope was found for the ternary system. 
All the ternary mixtures experimentally studied 


Cnao and O. A. 


fie and fi 


HovuGen 


boil at a higher temperature than the azeotropic 
mixture of cyclohexane and ethyl acetate. This 
binary mixture is indeed the lowest boiling of all 
the solutions prepared with the three given 
components. Fig. 11 that 
definite trend for all the equilibrium tie lines on 


shows there is a 
the ternary diagram to converge toward the com- 
position of this binary azeotrope. Data point 
ABC-79 the 
purpose of ascertaining if any ternary azeotrope 


was determined specifically for 
in the neighbourhood of this binary 
might form, but the finding was negative. The 


liquid phase of this ternary mixture contains only 


azeotrope 


2.9 mole per cent of benzene, but the boiling 
the 


equilibrium vapour contains even less benzene. 


point is thereby slightly increased, and 


THEORETICAL CONSIDERATIONS 


Equilibrium between the vapour and liquid of 
a multi-component system is attained when the 
pressure and temperature is the same throughout 
the system and when the fugacity of each com- 
ponent in the two phases is the same, that is 


fig =f. Under these conditions the relation of 


vapour composition to liquid composition for 
ideal solution behaviour is expressed as : 


(1) 


where at the temperature and pressure of the 


Yi Sic 


system, 


fugacities of pure component 7 in 
£ 
the respective vapour and liquid 


phases. 


Any lack of ideal solution behaviour in each 
phase is allowed for by arbitrary terms designated 
and y,,, thus : 


as activity coefficients, y,, 


Vix Vilic Yun 


Here the standard state fugacities f,, f,, are 


related to pressure by fugacity coeftlicients 
vj, and v,,, thus : 


(3) 


Sig = Vig and, 
total pressure 


P, = vapour pressure of pure liquid component # 
at the temperature of the system 
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At atmospheric pressure and below, the activity 
coefficient in the vapour phase y,, is unity and 
the ratio of fugacity coefficients v;_/v,, is nearly 


unity, hence 


ip 


Yip (4) 


and the 
activity coefficient y, is evaluated. 

In a binary system at constant pressure and 
temperature the 
efficients with composition is related by the 
Gibbs-Duhem equation : 


From measurements of y;, 2 


variation of activity co- 


[(2 In ¥,)/(d Ina,) = (d In y,)/(9 (5) 


The excess free energy »* of any component i 
in solution is related to its activity coefficient, 


— RT In y, (6) 


thus 


For a binary system : 


G® =a, + = RT (a, In y, + (7) 


For conditions of constant temperature and 
pressure Repuicu and Kistrer [13] have related 
excess free energy to composition by a series 
function using terms sufficient to fit the experi- 
mental data, thus : 


{G® = RT 2,2, [B’ + C’ (a, — 24) 
D’ (a, +.. bor (8) 


The terms 2, x, provide for the zero value of G® 
at compositions corresponding to the separate 
Formulations for individual 
activity coeflicients are obtained from equation 
(8) through the definition of a partial property, 
thus : 


pure components. 


[(d (ry Ng) G* RT y, (9) 
to give : 

2, 

(a, 

a, [B' (a, — 2,) 
D’ (a, - 
a,2,[B + 

D’ (a, 

a, [B’ (2, — 2,) + C’ (6x, 2, — 1) 
D’ (a, — (8,2, 1) +...] 


In yy 


C’ (6x, x, — 1) 
1)+...] 


(BX, 


Subtracting the two parts of equation (10) 

gives : 
In (y,/y¥e) = (v2 — 
+ D' (a —2,)(1 


(67, — 1) 


82, 4 (11) 


The same relationship for In (y,/yg) may be 
obtained by differentiation of equation 8 with 
respect to 


(G*/RT) 


In yy 
Y1 ve) lpr (12) 
Thus experimental data on activity coefficients 
can be correlated by using a single equation in 
terms of the ratio of activity coefficients without 
requiring separate expressions for the individual 
coefficients. 


IsoBparRiIc CONDITIONS 


For vapour-—liquid equilibrium in a_ binary 
system at constant pressure, temperature is fixed 
by composition and is not an independent 
variable. Hence the functional form of equation 
(8) may also be used to express excess free energy 
under isobaric conditions using slightly different 
The the 
equation are then dependent on pressure only. 


numerical constants. constants in 


Thus, at constant pressure 
{G® = RT 2, 2, [B 
D(a, —%)+.. 


the Gibbs—Duhem 


conditions of 


(13) 


In its conventional form 


equation is constant 
temperature and pressure. IsBL and DopceE [6] 
have extended the Gibbs—Duhem equation to 


and 


given at 


conditions of constant pressure variable 
temperature. Since In y is at constant pressure, a 
continuous and single-valued function of tem- 
perature and composition, for a binary system 
in the liquid phase, the following exact differential 


equations may be written :— 
[dln y, dT) ,., 


(d In Y1 pr 
(d In y,/d Ina,),7 4 In ay], 


(14) 


and similarly for component 2. 
The variations of activity coefficients with 
temperatures at and 


position are given as follows : 


constant pressure com- 
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(H, H,) (15) 


(d In V1 T) 


pri 


and similarly for component 2. Combining 
equations (5), (14), and (15) gives : 


[(d In y,) (dIna,) = (d In yg) (d In Z|, (16) 


where 
Z — (AH RT*) (dT dz,), 
AH =H —7,H, 
H, H,, H, = the molal enthalpies of mixtures 


and components, respectively, 
each in the liquid phase. 


The term AH is the integral heat of solution in 
dissolving liquid components in their standard 
states to form a liquid solution of given composi- 
tion at constant pressure and corresponding 
temperature. 

From equation (16) 


¥, +a diny, = Zda,), (17) 
From equations (7) and (16) : 


[(d /dx,) (G®/RT) 


In y, — In ye + (d In y,)/(d Ina,) (18) 
(d In y,)/(d In 


In y, — In + 


The values of Z are diflicult to determine experi- 
mentally since they involve measurements of the 
variation of temperature with compositions at 
constant pressure and the heats of solution 
which involve hypothetical standard states. These 
measurements are circumvented by including the 
effect of Z in the empirical equation for the ratio 
of activity coefficients, thus retaining the same 
form as for isothermal conditions but with con- 
stants having slightly different numerical 
values depending upon the magnitude of Z, thus ; 
In = O(a, + (62,7, 1) 

d (a, — z,)(1 — 87,27.) +... (19) 


Using identical forms for isothermal and iso- 
baric conditions is not mandatory. This proved 
satisfactory for the three binary systems under 
investigation, but it might require additional 
terms for systems where Z is large. With this 
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particular form, the integral In (y,/yq) is 


zero for both isothermal and isobaric conditions. 
For isobaric conditions, the integral is not 
necessarily zero, and might be accounted for by 
additional terms dissimilar to those in the iso- 
thermal equation. 

Subtracting equation (18) from (7) : 


In yg = (G®/RT) — a, [(d/dx,)(G®/RT) — (20) 


Combining equations (13), (19) and (20) :— 


In yg = 2, + C(x, — 24) 


D(a, z,)° 
(21) 
x, (Ua, — + (62,2, — 1) 
d (a, Ty) (Rr, 1) 
Similarly, 
In y, =a, 2, [B + C(x, — 
D(a, — 
+ [b (x, + (62, 1) 
d(x, Tq) (Br, 


In calculating vapour compositions, tem- 
perature sensitivity is minimized by using ratios 
of activity coefficients and vapour pressures by 
combining equation (4) for the two components, 
thus : 

= (71% (ye Po) 


= 1/E (22) 
E =(y, Ps) + 1 


For calculating temperatures under isobaric 
conditions it is necessary to have equations for 
the separate activity coeflicients, and then by 
trial and iteration to establish the temperature 
which satisfies the relation y,7 =a, y, Py; in 
agreement with the correct value of the vapour 
pressure of the pure component. 


TERNARY SYSTEMS UNDER ISOBARIC 
CONDITIONS 
The excess free energy of a ternary mixture ts 
related to the activity coefficient by the expres- 
sion : 


gE E E 
Gizg = 214) + + 


= RT (a, In y, + In yg + agin ys) (23) 


_ 
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Under conditions of constant pressure : 
G3) 


RT x, (B 


Gs 
— 


Here the first three terms on the right side 


(24) 
C2 


represent contributions by the separate binaries 


ry 
The 


assigning 


and the last term represents ternary effects. 


binary contributions are related by 
aggpuapente constants to equation (13), thus for 


G5: 
RT x, | By 


Cy2 - 


(2 


+...) (25) 


'E 
Gio 


As the temperature is fixed by composition the 
constants in equations (24) and (25) are dependent 
on pressure only. 


The 


ternary system under isobaric conditions is given 


modified Gibbs-Duhem equation for a 


by an expansion of equation (17) thus : 


a, diny, + din yg + In 


(AH/RT®) dT (26) 


On differentiating equation (23) and combining 
with (26) the individual activity coefficients may 
be expressed in terms of G® 


In yy 


AH {dT 
an (22) 


7+ 

a, iF (a7) 
{GE 

a, 


for the other components. 


021/29 


The 


and similarly 
ratio of the activity coeflicients follows :— 


In (71 /¥2 
[(d/d (G®/RT)) 
(AH /RT*)(dT/d 


r (28) 


and similarly for the other ratios. 
Expressions of activity coefficients and their 
ratios explicitly in terms of compositions result 


upon substituting equation (24) into e quations (27) 
The difficulties with the 
(AH /RT*) (dT 
modifying the constants as done 
with the binary system, thus : 


and (28), respectively. 


correction terms are cir- 


cumvented by 


-(G,E,/RT) + (a +2 ) 
[bos %)+...] 

[bs, C31 (V3 

2, (7, + (Cog + 
[ 272, (7, + New +...) 

(2, ls 


(0 
Cy 


(2c, 


In y, 


C23 (12 


Ys) 


+ 


where G,£, is given by equation (24). 


In (y,/¥2) 


(30) 


Similarly expressions for the other components 


are obtained by cyclic advancement of the sub- 
scripts in the order 1, 2, 3, 1. 

In calculating vapour composition from these 
equations, temperature sensitivity is minimized 
by using ratios of activity coefficients and vapour 
pressures by combining equation (4) for the 


three components, thus :— 
y, =1/E 
Ys = Ps/¥, 7, P,) E 
1 


E - r 2 P,) (47 P,)} 
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0-4 


Mole fraction of 


Fic. 7. 


Binary DaTa AND CORRELATION 


A total of 13 
each of the three binary systems at 760 mm Hg 


t values were measured for 


pressure and for nearly equal intervals of liquid 
composition. These are presented in Tables 5, 6, 
and 7 and Figs. 7, 8, 9 and 10. For each binary 
system the data were correlated with a series of 


° 


Ethyl acetate (1 
~ Benzene (2), 


Benzene (1) Cyclohexane (1) 
~Cycionexane (2)—9% -Ethyi acetate (2) 
4 
A 


° 


° 


Mole fraction of component 1 in vapour 


+O 


02 O68 O8 
Mole fraction of component 1 in liquid 


Fie, 8. Equilibrium compositions in the binary systems. 


Cuao and O. 


Benzene (!)-Cyciohexane (2 


A. 


HouGEN 


2) 


acetate (1) -€ 


component 


Equilibrium temperatures and compositions in the binary systems. 


Redlich-Kister of 


plexity, namely, with 1, 2 and 3 constants. 


com- 
These 
the 


method of least squares using all experimental 


equations increasing 


constants were in each case evaluated by 


data. 


The 


equations 


following constants were evaluated for 


relating activity coeflicient ratios, 


using logarithms to the base 10 : 


Ethy! 
acetate (2) 


Ethyl acetate (1)- Benzene(2) 
03 04 O6 OF O8 OD 
Mole fraction of component | 


JS Benzene (1)- Cyclonexane (2) 


Fic. 9. Excess free energy of the binary mixtures. 
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Mole fractior 


Ratio of coeflicients in the binary 


activity 


systems. 


Table 5. 


2-Constant | 3-Constant 


1-Constant 


For ethyl acetate [1 
benzene [2 ‘ 0-0041 
d 00-0028 


0-042 


For 60-1454 O-l444 
cyclohexane [2 O-O104 
d 0-002 


benzene [1 O-1445 


For cyclohexane [1] b O-3B850 O-3849 O-3848 


ethyl acetate [2] 


d O-O076 


Calculated values of vapour composition were 
compared with experimental values for each 
binary system, using each correlating equation. 
The root mean square deviations of calculated 
values from experimental values of vapour com- 
positions expressed in mole fractions were 
obtained (Table 6). In these calculations vapour 
pressure ratios were based upon experimental 


equilibrium temperatures. 


Table 6. Root mean square deviations of 


(Y cale. Yexp.) 


1-Constant | 2-Constant | 3-Constant 


equation equation equation 
Kthvyl acetate 
benzene 0-0007 


00-0022 


O-001L0 


cyclohexane 0-012 0-0010 


ethyl 


Benzene 
Cyclohexane 


acetate 00-0059 00-0020 


the above deviations 


indicates that for the ratio of activity coefficients 


An examination of 


an optimum correlation is obtained with 2- 
constant equations for the benzene — cyclohexane 
and cyclohexane —ethyl acetate systems, and 
that a 1 


acetate 


constant equation suffices for the ethy! 


benzene system. By optimum cor- 
relation is meant that the use of additional terms 
does not give significant improvement to justify 
the extra work involved. 

For simultaneous correlation of vapour — liquid 
composition with equilibrium temperature at 
constant pressure a trial and iteration procedure 
was employed. Calculation of separate activity 
coeflicients as well as their ratios were required 
for this purpose. The procedure followed was to 
calculate vapour pressures of the separate pure 
components at an assumed temperature from the 
The 


position calculated thus from equation (22) 


vapour pressure equations. vapour com- 
was 
subsituted into equation (4) to verify the assumed 
vapour pressure ; the second temperature to be 
assumed was then obtained from this vapour 
pressure through an equation for vapour pressure. 
The process was repeated until successive trials 
gave the same temperature within 0-01°C, 

The following equations were obtained for 
individual coeflicients in 


expressing activity 


terms of logarithms to the base 10 :— 
#, [B 

a, [B 


a, [b (a, — a) 


log y C (a, 


C (a, 


log Ye 


A comparison of calculated with experimental 


££. 
\. Cyclohexane (1) -Ethy! acetate 
O-2 
O-1 Benzene (1) 
| 4 
as 
\ 
| — 
\ 
VOL, - \ 
957/58 | 
O 0-2 0-4 
| 
— 
1)] 
(33) 
7, — 1)] 
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For ethyl acetate 


[1] — benzene [2] | 0-0372 0 0-0359 0 
For benzene [1] - 
cyclohexane [2] 
For cyclohexane 
[1] — ethyl ace- 


tate [2] | 0-3849 | 0-0321 | 


| 
| 
| 


0-1443 | 0-0100 | 0-1507 | 0-0100 
| | 

| 


0-4190 0-0149 


values for these three equations is shown in Figs. 
7, 8, 9 and 10. 

The root mean square deviations of calculated 
from experimental values of temperature and 
vapour composition are tabulated herewith. A 
comparison with deviation obtained from the 
conventional method of correlation is also given. 
The conventional method does not recognize the 
the the 
and activity 


constants in free 
the 


coeflicient ratio expression. Those constants are 


difference between 


energy expression those in 


Table 8. 


Veale Yexp. 
mole fractions 
Con- Con- 


System Isobaric | ventional | Isobaric | ventional 


Ethyl 
benzene 
Benzene 


acetate - 


00-0007 0-0007 0-07 0-07 


cyclohexane 0-0012 0-0012 0-06 0-13 


Cyclohexane - 


ethyl acetate 0-0022 | 0-0022 0-13 0-37 


then solely determined from the expression of the 
activity coefficient ratios as in equation (19). 

It will be seen that the correlations are nearly 
the same with isobaric and conventional equations, 
The calculated temperatures deviate from the 
experimental values by about 0.1°C for all three 
systems. Employing two sets of constants did 
not improve the calculated values of vapour 
composition but the cal- 
culated values of temperature. 


slightly improved 


Cuao and O. 
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TERNARY SYSTEMS 
Experimental values of vapour-liquid com- 
positions at atmospheric pressure are shown by 
the lines connecting the vapour with the equili- 
Values of 


x, y, and ¢t were measured at 760 mm Hg pressure 


brium liquid composition (Fig. 11). 


for 77 ternary compositions uniformly distributed 
over the entire range of ternary compositions. 
These were in addition to the 39 compositions for 
the three binary systems, and were recorded in 
Table 14. 


Table 9. Constants for the ternary system for use 


in expressing log y, ‘ye 


3-Constant | 6-Constant | 9-Constant | 


equation equation 5-Constant 


(1 for each | (2 for cach 
pair) 


equation 
(3 for each | equation 


pair) pat) 


| ©0372 
00-1454 
0-3850 


00-0369 
O- 1444 
0-3848 
0-0042 
0-03518 


00-0372 
00-3849 


0-0369 | 
0-1443 
0-3849 | 
0-0041 
0-0010 
0-0321 


CAB 
boa 
“BC 


0-0100 
00-0321 


00-0029 
0-0076 


dnc 
dea 


dap | — 00028 


Correlations of ternary data were made similarly 
to those for the three binary systems, basing them 
first upon experimental measurements of tem- 
perature and then correcting the temperature to 
give complete correlation for all equilibrium 
conditions. The first correlation was made with 
experimental temperature values, using 1-, 2- and 
3-constants for each binary pair with corres- 
ponding 3-, 6- and 9-constants for the ternary 
solution. No additional constants were required 
for ternary effects. 

The root mean square deviations of y.... 
Yexp, in mole per cent are shown in Table 10. 

The accepted equation for relating all three 
activity coeflicient ratios in the ternary system 
5-constants; for individual 


required only 


activity coefficients at constant pressure, 5 


permitted accurate cal- 
temperature. With 


equations (29) and (30) for activity coefficient 


additional constants 


culation of reference to 
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Table 10. 


Table 12. 


| 9-Constant 5-Constant 


3-Constant 6-Constant 


equation equation equation equation 


0-13 


0-13 0-16 


ratios and activity coefficients, the following 
constants apply. They are given for use with 
logarithms to the base 10 instead of natural 
logarithms. 

Table 11. 


| 


bie 0-0372 By 0-0359 
bes 01443 Bos 0-1507 
bs, 0-3849 Bs) 0-4190 
0 Cis 0 

Ces 0-0100 Cog 0-0100 
00821 Ca 0-0149 


designated as 


Here 1 


component 


ethyl 
A), 2 


3 = cyclohexane (component C). 


acetate (also 


benzene (component B), and 


The root mean square deviation of calculated 


values from experimental measurements are 
presented in Table 12. 

It will be seen that the correlations are nearly 
the same with isobaric and conventional equations. 
The temperatures calculated with the isobaric 


equations deviate from the experimental values 


| Veale. ~ Yexp. 
(mole 
| per cent) 


0-13 
0-13 


Isobaric equation 
Conventional equation 


by about 0-1°C, Employing two sets of constants 
did not improve the calculated values of vapour 
composition but slightly improved the calculated 
values of temperature. 
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The error in thermodynamic predictions as a function of the 
displacement from equilibrium 


R. P. Rastoci and K. G. Denpicu 


Dept. of Chemical Technology, West Mains Road, Edinborough 9 
(Received 17 September 1957) 


Abstract—The paper is concerned with the extent to which a thermodynamic prediction of 


the effect of temperature or pressure on the yield of a reaction is in error, due to the non-attainment 
of equilibrium. The magnitude of the error is evaluated as a function of the fractional attainment 


of equilibrium and is shown to be surprisingly large under conditions such as might be encountered 


in industrial practice. 


Résumé 


Cet article traite des écarts intervenant dans les prédictions thermodynamiques des 


effets de température ou de pression sur le rendement d'une réaction quand l’équilibre n’est pas 
atteint. L’auteur évalue lerreur en fonction de approche de l’équilibre et il montre que dans 
le cas des installations industrielles elle est étonnament grande. 


Zusammenfassung 


Voraussage des Temperatur- oder Druckeinflusses auf die 
sein kann, bedingt durch die Nichterrei hung des Gleichgewichtes. 


Die Arbeit beschiiftigt sich mit der Frage, wieweit eine thermodynamische 


Ausbeute einer Reaktion fehlerhaft 
Die Grdésse dieses Fehlers 


wird als Funktion der Annd&herung an das Gleichgewicht errechnet und erweist sich als fiber- 


raschend gross unter Bedingungen, wie sie 


One of the useful applications of thermodynamics 
is the prediction of the effect of changes of 
temperature, pressure or concentration on the 
The 


equations used for this purpose are based on the 


state of physical and chemical systems. 


second law and are thus strictly applicable only 


if there is true equilibrium—a condition never 
actually attained in industry, due to the time 
What is the error 


in these predictions, as a function of the degree 


factor. The question arises : 
of attainment of the equilibrium state ? 

This problem, which it seems has not been 
discussed previously, is obviously of greater 
significance in chemical than in physical processes. 
In ammonia synthesis or sulphur dioxide oxida- 
tion, for example, the chemical engineer uses Van’t 


Hoff’s equation, 


dink AH 


1) 
dT RT? 


for the purpose of predicting the effect of tempera- 
ture on the yield of the reaction. This prediction 
is inevitably somewhat in error whenever the 
contact time in the plant is less than infinite. A 
more exact prediction must be based on chemical 


in der industriellen Praxis auftreten kénnen. 
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kinetics, which allows for the finite speed of 
the forward and backward reactions. This can 
be done only when the form of the kineties is 
known, together with experimental data on the 
temperature coeflicients of the velocity constants. 

In the shall the 


prediction based on thermodynamics with the 


present note we compare 
more exact prediction based on chemical kinetics. 
For this 
reactions whose kinetics are sufficiently simple. 


purpose it is convenient to examine 


Consider the reversible reaction 
A+B X+Y 


assumed to be second order in 
The reagents A and B will be 


taken as being present in equimolar amounts, 


which will be 


both directions. 
their initial concentrations being a,. Let x and 
r, be the concentrations of X (or Y) at time ¢ 
and at equilibrium respectively. The reaction 
vield at time ¢ is 


¢=21/a, (2) 
and the equilibrium yield is 


>, = 


‘ 
VOL. 
= 
(3) 
|| 


R. 


P 


Assuming that the system forms a_ perfect 


gaseous mixture, the equilibrium constant is 


2 
r 


(4) 


(do 


(5) 


Hence by differentiation and using equation (1): 


dd, AH 
dT 2RT?(i + Kip 


and this equation shows how the yield varies 
with temperature if there is true equilibrium. 
Consider the matter now according to chemical 

kinetics : 

da 
k, (ag — kya (7) 
where k, and k, are the velocity constants of the 
forward and reverse reactions respectively. In- 


tegrating we obtain 


t 
2h 
in | x z x x (8) 
\\Ki + \Av —a — Ka, + 


Hence solving equation (8) for the ratio ¢ = 2, ay, 
which is the yield at time t, we obtain 


exp (2k, a,t K 


kt + 1) exp (2h, ao t 


(9) 


Differentiating with respect to temperature and 
using equation (1), together with 


dink, &E, (10) 
dT RT? 
where E, is the activation energy for the forward 
reaction, we obtain, 
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(= 4a, k, t(2E, — AH) exp (2k, K+) — K! AH [1 — exp (4k, ay t K~4)] 
2RT? —1) — + 1) exp (2k, ag t 
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From this equation the time ¢ can be eliminated 
by substitution from equation (8). Let the 


fractional attainment of equilibrium be defined as 


f (12) 
Then we finally obtain 
dd 
(37), 


le 


+ (2E, — AH)(1 — f)(1 +f — 2f4,) In 


+f 
1—f 


is related to AK through equation (5). 


(13) 


where ¢, 
This equation shows how the yield varies with 
temperature as a function of f, and it reduces to 
equation (6) when f = 1 — Le. at thermodynamic 
equilibrium. 

The ratio r of the true temperature coefficient 
of the yield to the temperature coefficient which 
would be predicted on the assumption that the 
system is at equilibrium, is thus obtained by 
dividing equation (13) by equation (6) and is 


_ (d¢/2 T), 
dd, dT 


‘(1 (1 
(1 \J ( ( f 


> (14) 
[1 — 24,f(1 — f) — In 
AH * 


(24. — 1)]) 
1—f 


This equation is also applicable to reactions 


whose stoichiometry is 

provided that these also have second order 
kinetics. A well-known example for which 
experimental data are available is the decomposi- 
tion of hydrogen iodide : 


2HI = H, +I,. 


(11) 


4 
or 
(1 — ¢)? = 
1957 
r 
Zz 


The error in thermodynamic predictions as a function of the displacement from equilibrium 


For this reaction the equilibrium constant, K, 
has a value of 2-18 « 10° at 763-°8°K* and E, 
and AH of 44,000 cal/mole and 
3000 cal/mole respectivelyt. 


have values 


Fic. 1. The ratio r as a function of the fractional attain- 
ment of equilibrium in the reaction 2 HI Hy + Ig. 


These figures have been used in a numerical 
solution of equation (14). The results are shown 
in Fig. 1, where the ratio r is plotted as a function 
of f, the fractional attainment of equilibrium. 
It will be seen that, for a value of f of 0-9, the true 
temperature coefficient of the yield is 10-6 times 
the value predicted thermodynamically on the 
supposition of equilibrium. When the system is 
1 per cent from equilibrium (f = 0-99) the corres- 
ponding factor is still as much as 2-8. In this 
instance, at least, the effect of a few per cent 
displacement from equilibrium clearly gives rise 
to a disproportionately large error in the ther- 
modynamic prediction, The for this is 
that the temperature coefficient ‘of the equili- 
brium yield is determined by AH, whilst the 
temperature coeflicient of the yield when the 


reason 


reaction is far from equilibrium is mainly deter- 
mined by £,, the activation energy of the forward 
reaction. Thus whenever E, and AH differ by an 
order of magnitude, as in the case of the decom- 
position of hydrogen iodide, the ratio r may be 


expected to be fairly large at comparatively 


small departures from equilibrium. 

This has been confirmed by solving equation 
(14) for various values of the ratio EF, AH at a 
fixed value of f equal to 0-9, and for fixed values 


*Tayvor A. H. and Crist R. H. J. Amer. Chem. Soc. 
1941 63 1376. 

+For a reaction of the type 24 = X 
ship between K and ¢, is 4K = ¢,7/(1 
equation (5). But equation (14) applies unchanged. 


Y, the relation- 
¢,)*, in place of 


3 


Influence of the relative magnitudes of E and AH. 


of 
expected, the ratio r has a large (positive or 
E, and AH differ 


Higher values of ¢, also tend to 


The results are shown in Fig. 2. As 


negative) value whenever 
considerably. 
increase the magnitude of r. 

In addition to the temperature coefficient, 
dé/dT, it is of the 


pressure coeflicient, dé dP, and the extent to 


interest to consider also 
which its equilibrium value differs from its value 
prior to the attainment of equilibrium. — For 


this purpose consider a reaction of the type 
A 


Fx wr 


will be supposed that the reaction 


taking place in a perfect gas mixture. 
simplicity, it 
takes place in a closed system of constant volume 
(and temperature) rather than in a flow system. 
For this reason the pressure P varies during the 
course of the reaction and is not, properly speak- 
Instead we shall 
the 
or, which comes to 


ing, an independent variable. 
P,, the 


before reaction commences 


choose initial pressure in vessel 
the same thing, the initial concentration a, of 
the two reagents, taken as equal. 

Let é and ¢, be the yield of product at time 
t and at equilibrium respectively. The variation 
of the equilibrium yield with the initial concentra- 
tion is given by 


dd, 


(1 
Ay (1 T 


day 
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be=0'5 
| 
| 0-025 | 
4 | | 
| | | | 
4 | | 
r 
| 
6.2, 
VOL, 
957/58 
| (15) 


R. 


where ¢, is calculable from the equilibrium 


constant : 


(16) 
a, (1 — ¢,)* 

If the reaction kinetics happen to correspond to 
the stoichiometry, we have 


dd 2 


ke ® (17 ) 
dt 


which may be integrated to give ¢ as a function 
of t. Proceeding as in the preceding example, 
we finally obtain r, the ratio of the true dependence 
of the yi ld on the initial concentration, to the 


thermodynamic value as given by equation (15): 


dy), 
dd, da, 
., ) 
f (18) 
(1 — 
(1—fé,”) 
In | | fb, (1—¢,°) 


where f, the fractional attainment of equilibrium 


is defined as in equation (12). 


Fic. 3. The value of r in a reaction influenced by the 
pressure. 


Using equation (18), the value of r has been 
plotted in Fig. 3 as a function of ¢, (which is 


determined by A and a,), for a fixed value of f 


equal to 0-9— which is about as high a value as 


and K. G. Densicu 


might be expected under normal industrial 
conditions. It will be seen that as ¢,, the equili- 
brium yield, becomes greater the larger the error 
in the thermodynamic prediction. This is as 
might be expected : in a reaction of stoichiometry 
A+B X the effect of an increase in pressure, 
or concentration, is to raise the rate of the 
forward reaction more than that of the reverse 
reaction, and thus to increase ¢. If ¢, is large, 
the effect is almost entirely due to the forward 
reaction, the backward reaction being relatively 
insignificant. Thus, under these conditions, a 
prediction based on equation (15), which implies 
that the forward and backward processes are 
proceeding equally, would be a great deal in 
crror. 

The result of this short survey is to show that a 
thermodynamic prediction can be considerably 
in error when a chemical system Is, say, 10 per cent 
away from equilibrium (f = 0-9). In the case of 
the temperature coeflicient of the yield, especially, 
the error can b quite large whenever the ratio 
Ail Is large. 

The comparison has been extended to other 
types of reaction, of different kinetic orders, and 
similar results have been obtained. The general 
conclusion to be drawn is that thermodynamic 
predictions are to be used with a good deal of 
caution, whenever the system in question is 
appreciably displaced from equilibrium. Of 
course the great advantage in the use of thermo- 
dynamics is that it requires a minimum of 
experimental data—and also it shows the 
limiting value of what is attainable. To apply 
kinetics requires a much larger amount of 
information, in particular a knowledge of the 
kinetic order of the reaction and also—for a 
prediction of the effect of temperature—a 
knowledge of the activation energies. Nothing 
that has been said in this short paper need be 
taken as diminishing the great value of thermo- 


dynamics in chemical engineering. 
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A note concerning turbulent exchange of heat or mass with a boundary 
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(Received 18 September 1957) 


Abstract 
of the equation 


A renewal model is proposed which allows the derivation, at large Prandtl! numbers, 


Nu ~ Re Pri, 
2 


where Nu, Re and Pr are Nusselt’s, Reynolds’ and Prandtl’s numbers, and / is the friction factor. 


Resume 
léquation : 


Lauteur propose, dans le cas des nombres de Prandtl élevés, une forme nouvelle de 


Nu ~ Re Pr} 
Aj 2 


propre 4 la dérivation ot, Nu, Re, Pr, sont les nombres de Nusselt, Reynold et Prandtl, et f le 


coellicient de frottement. 


Zusammenfassung—Es wird ein Austauschmodell vorgeschlagen, das die Erweiterung der 


Gleichung 


Nu ~ Re Jf Pri, 


auf grosse Prandtl-Zahlen gestattet. Hierin bedeuten Nu, Re und Pr die Nusselt-, Reynolds- 
und Prandtl-Zahlen und den Reibungsfaktor. 


IN their recent papers Hanrarry [1] and Jounson 
and Hvuane [2] extended the model of continuous 
renewal of the interface layer, as given by 
Danckwexts [3], to the case of mass exchange 
with a boundary. A similar idea has been used 
independently by the author in [4], with a more 
complete result than the one obtained in [1] 
and [2]. 

In these papers it is assumed that, owing to 
turbulence, the liquid layer near the boundary is 
continuously replaced by other liquid layers, 
each of them remaining at the boundary for a 
very short time. 

Observations made by Face and TowNnenp [5] 
in the immediate vicinity of the wall on a turbulent 
fluid in which they traced—-by the aid of an 
ultra-microscope—the motion of dust particles, 
perhaps provide experimental support for a 


renewal mechanism. These authors noticed the 
existence of large lateral displacements of groups 
of these particles, and mentioned the fact that 
the motion can be regarded as practically rec- 
tilinear only during the interval between two 
such displacements. This experimental evidence 
is in contradiction with the model of the laminar 
sublayer which implies the existence of a region 
within which the motion should be rectilinear, 
but is consistent with a renewal mechanism. 

In [4] the author obtained for the coefficient 
of heat transfer at large Prandtl numbers, in 
the case of a tube, 


Nu ~ Re Jf Pr, (1) 


Nu Nusselt’s number, Re Reynolds’ 
number, Pr Prandtl’s number and f the friction 


being 
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factor. In deriving equation (1), dimensional 
considerations were used besides DaNcKWERTS’ 
renewal mechanism. 

Within the framework of classical theory one 
obtains for Nusselt’s number, at large Prandtl 


numbers [6], [7], 


Nu ~ Re J prs (2) 


It should be stressed though that (2) has been 
derived on the assumption that in the close 
neighbourhood of the wall the eddy diffusivity 
is proportional to the third power of the distance 


y | 
~~ Vp (3) 
\ 


Here v is the kinematic viscosity, y the distance 
from the wall, 7, the shear stress at the wall 


from the wall : 


and p the density. This assumption is needed to 
obtain agreement with experimental data. 

In this note a renewal model differing in some 
respects from the one used in [1], [2], and [4] 
is proposed ; thus equation (2) may be derived 
without making use of the above assumption. 
It will be assumed that the boundary is divided 
along the direction of flow into short sections of 
length a. A certain part of the liquid flows 
along the wall a distance a), after which it mixes 
with the bulk of the liquid. This process 1s 
repeated at intervals 2, apart. Let it be assumed 
that the liquid layer which flows along the wall 
in each section is non-turbulent and that the 
* path-length ” 2, is small enough to render the 
thickness of the laminar boundary layer corres- 
ponding to the flow along a plate smaller than the 
thickness z, of the liquid layer. Then the velocity 
distribution in this layer can be approximated 
by the velocity distribution in the case of a flow 
along a plate [8]: 


2.4 
Uo a 7 
1-33 2\4 
— (“2 = 1-78 

4 vz 
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Here u, is the velocity of the liquid layer at the 
beginning of each section, x is, for each section, 
the distance from its beginning, and 6 the 
thickness of the boundary layer, given by : 

5 = 52 (— (5) 


Ny 


The equation of heat transfer in a liquid layer 
which flows along the wall in a section ts : 
u + (6) 
oY oy* 
where T is the temperature and a the thermal 
diffusivity. 

If Prandtl’s number is large enough, the tem- 
perature does not vary sensibly except close to the 
wall, and we can assume the temperature of the 
liquid layer which flows at the wall to be the 
same, for z 0, to the average temperature. If 
the thickness oy of the thermal boundary layer 
(under conditions of laminar flow) is smaller 
than the thickness of the layer flowing along the 
wall in a section, the temperature distribution 
may be approximated by one for a semi-infinite 
medium. (Since 6, < 8 if Pr > 1, it is evident 
that for Pr l o7 z, if also 6 z,). Hence, 
we have the following boundary conditions : 


T=T, fory=0 
T=T,, fory= (7) 
T = T,, for z = 0. 


The solution of equation (6) with boundary 
conditions (7) has been given by PoHLHAUsEN 


[8], [9], [10] : 


T,, —T, 
( exp { — 0-22 Pr 2°} dz) / 


0 


[exp {— 0-22 Prz*}dz. (8) 
Since 


a 


0-89 
exp { — 0-22 Pr zy} dz 


(0-22 Pr)’ 


| 
VO 
1957 
uy 
4: 
(4) 
2\4 
v = 4b 
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equation (8) becomes 
= 
(0-22 


‘xp { — 0-22 Pr 2*} dz. 9 
O-89 { (9) 


For the heat transfer coefficient we obtain 
easily : 


h, = 0-34 Pr. (10) 
va, 
The average value h of h, is given by : 
= hyde. (11) 
Zo 
Hence, 
0-68 ) Pri (12) 
v To 


Taking into account that the “ path-length ” 
2 is very small, h may be regarded as a quasi- 
local quantity, the quasi-local coefficient of heat 
transfer. 

The “ path-length ” a2, as well as the velocity 
Uy, are In other 
words 2, and are functions only of p, and 
which characterize the hydrodynamical state 
of the fluid in the vicinity of the boundary, and 


hydrodynamical quantities. 


not of A, for instance. Dimensional considerations 
permit writing : 


~ (13) 
~ (14) 


Taking into account (13) and (14), equation 
(12) becomes 


(15) 


Using the expression, 


where u,, is the average velocity, and the usual 
dimensionless groups, equation (15) becomes : 


Nu ~ Re J 4 Pri, (16) 
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In what follows the assumption that z, is 
constant in time will be abandoned ; for a given 
value of x), equation (12) will be used to give the 
approximate coefficient of heat transfer. Con- 
sider an infinitesimal region chosen at random at 
the boundary, and let ¢dz, be the number of 
cases in which this region belongs to a path- 
length between 2, and a, + dry. 

Dimensional considerations lead to : 


a 
where 
(18) 


It can easily be shown that 2, is proportional 
to the average value of 2p. 

The quasi-local coefficient of heat transfer h 
can be defined as the average value of h. 

This treatment of the problem holds true only 
if ¢ takes finite values in a sufficiently narrow 
interval about the average value of 2, and 
practically vanishes outside this interval. Only 
in this case can the motion in the vicinity of 
the wall be considered as quasi-stationary, and 
hence equation (12) may be used. Using equations 
(12) and (17) we obtain 


-4 
h 0-68 A (19) 
| F a2 
a) 


Since the ratio of the integrals is a constant, 
by using equations (13), (18) and (19), an equation 
of the same form as equation (16) is obtained 


for h. 


NOTATION 


a = thermal diffusivity 

jf = friction factor 

T = temperature 
T, temperature at the wall 
T,, = average temperature 

u =2 component of velocity 
v = y component of velocity 
az = co-ordinate in flow direction 
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= distance from wall density 
coefficient of heat transfer = shear stress at the wall 
= eddy diffusivity Re = Reynolds number 
thermal conductivity Nu = Nusselt number 
kinematic viscosity Pr = Prandtl number 


REFERENCES 
Hanratry J. T. Amer. Inst. Chem, Engrs. J. 1956 2 359. 
Jounson A. I. and Huane C. J. Amer, Inst, Chem. Engrs. J. 1956 2 412. 
DANCKWERTs P. V. Industr. Engng. Chem. 1951 43 1460. 
E. Rev. Chim. 1957 in press (in Roumaniun). 
Face A. and Townenp H.C H Proc. Roy. Soc. 1932 A 135, 656 (quoted in [1)}). 
RUCKENSTEIN E. Stud. Cercet. Fiz. 1954 4 123 (in Roumanian). 
Risaup G. C. R. Acad. Sci Paris. 1956 242 959 
Levict V. G. Fiziko-khimicheskaia gidrodinamica. lad. Alead. Nauk. SSSR, Moscow, 1952. 
Broun E. A. Introduction a l'étude de la couche limite. Gauthier-Villars, Paris, 1955. 


POHLHAUSEN E, Z ingew. Math. Mech. 1921 1 115. 


h 
(1) 
[2] 
[3] 
(4) 
[5] 
(6) 
[7] 
[8] 
~~ | (9) : 
[10) 
vc 
1957 
268 


Discussion of the paper * Fugacities in simple gas mixtures "’ 


J. M. Prausnirz* 


(Received 3 May 1957) 


Tue paper mentioned above presents an interesting and 
useful method for estimating the fugacities of individual 
components in a particular type of gas mixture, and it is 
correctly represented as such, 

To illustrate the method, Prof. Prausnirz chooses the 
T data 
are available at 50°C [2]. He compares the results of 


Ne ( eH, system, for which extensive p t 


fugacity calculations by his simplified method with the 
smoothed results of Bennerr [1]. It is not pointed out, 
however, that the method used by Bennett for smoothing 
his results is based on the same general assumption of the 
Hence the 
comparison made is between two approximate sets of 
results. It 


type of mixture as is made by Prausnrtz 
is true, of course, that Prausnirz’ method is 
based on whereas 
Prof. method actual 
iD t T data. But both methods are limited by the 
assumption that the volume change of mixing is a sym- 
\ better test of 
Prausnitz’ method would be a comparison of his results 


a principle of corresponding states, 


makes use of the 


metrical function of the composition, 
with the more accurate values of activity coefficients for 
this system given in a paper by the writer [4 

Prof 
entirely 


were carried out by 
and the 
inaccurate. He 


Bennetr’s calculations [1 
differentiations 
fully recognized 
this difficulty, and pointed out that his calculated results 
did not meet the test 


made according to the 


graphical methods, 


involved are inherently 

of thermodynamic consistency as 
method of 

13 In order to get a consistent set of results from his 

calculated 

RepLicu 


Repuicu and Kisrer 


calculations, he smoothed his values by 


arbitrarily drawing straight lines on the and 


Kuster plot of In (yng VS. YNe 


Table 


Such straight lines correspond to a system of Type II 
as described by Repiicu and Kisrer [3], and are truly 
representative only of systems whose volume change of 
mixing curves are symmetrical functions of mole fraction. 
This is not strictly true for the Ng — Coll, system as may 
be seen from Bennett's Fig. 1. Prof. Bennerr clearly 
stated that his calculated values might be more accurate 
than his smoothed values even though thermodynamically 
inconsistent. The equations for In y which result when 
the straight-line relationship is assumed for the Repiicu 


and Kisrer plot are given by BENNETT as : 


In and In yg 


By? 


These are identical with the PrRausnitz’ equations (16) 
and (19). Prof. Prausnirz has presented us with a new 
method of calculating B (K in his paper). However, the 
comparison of his results with BENNeT?T’s smoothed values 
begs the issue of their accuracy. He is comparing two 
methods of calculating a constant for an equation which is 
not strictly valid for the system considered, although it is 
admittedly a good approximation over much of the range 
to which it was applied by Prausnrrz. 

The more accurate results calculated by the writer [4] 
test of 


smoothing, but did not 


met the thermodynamic consistency without 


result in straight lines on the 


Repicu and Kisrer plot. The curves that were obtained 
could certainly be approximated by straight lines over a 
considerable portion of the concentration range, and the 
method presented by Prausniiz apparently is equivalent 
to such an approximation in the pressure range to which he 
has applied it. However, there may well exist regions of 


concentration in which the straight-line approximation 
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does not give very good results. For example, in the 
No 
of both Bennerr and Prausnirz for yx. at low concen- 


Call, system at 100 atm the approximate values 


trations of nitrogen are appreciably different from the 
values determined by the writer. It should be noted that 
Benxnetr was unable to calculate the limiting values of 
y at infinite dilution of cither component. Hence his 
smoothed values of y at this condition r present extra- 
polated values. A comparison of y values for a pressure of 
10) atm obtained by the three different authors is 


tabulated below. 


Certainly the new method proposed by Prausnirz 
should be useful to those who require easy estimations of 
fugacities for the components of simple gas mixtures, but 
it should not be used with quite the confidence suggested 
by Prof. Prausnrrz’ comparison of results. 


Henverick C, Ness 


Department of Chemical Engineering, 
Rensselaer Polytechnic Institute, 
Troy, New York, U.S.A. 
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Chemical Engineering Practice. Edited by H. W. 
Editor Treror Davies). Vol. I: 
95s. Volll: Solid State vi 

Scientific Publications, London, 


Cremer (Managing 
General xiv 104 pp. 
632 pp., Butterworths 
1956. 95s. 
Vouume 
The Origins of Chemical Engineering 
The 


Economics of 


Chemical Engineer 

Production as Exemplified in 
Process Industries 

Investigation and Development of an Indust rial 
Process : 
Preparation of Materials Balances 
Preparation of Energy Balances 

Pilot Plants and Semi-Commercial Units 

The Preparation of Flow Diagrams for Full-Scale 
Production 

Units, Dimensions and Calculations. 


Vouume Il 


Fundamental Concepts of Matter in the Solid State 

Metals and Metallic Alloys 

Mechanical and Physical Properties of Plastics and 
Glasses 

Corrosion of Metals 

Porous Masses 

Porous Powder Metallurgical Products for Filtration 
and Related Uses 

Powder Metallurgy 

Important 

Transpiration Cooling 


Industrial Uses of Porous Masses 


Tue two volumes under review are the first of a set of 
twelve. The object of the series is stated in the Preface : 

“ Thus. commercial success in translating a laboratory 
method of a preparation into a full-scale manufacturing 
process depends as much upon the careful study of plant 
and plant design as upon consideration of the precise 
chemical reactions to be employed ; in short, industrial 
efficiency and the profits expected to accompany this can 
only be realized by sound chemical engineering. 

“« In order to achieve this end and to meet the intensive 
competition of the present age, the closest co-operation 
worker and the 


necessary 


between the more strictly theoretical 


technologist is essential. Furthermore, it is 
for both the pure and the applied scientist in industry 
to have an awareness of the significance of the economic 
engaged in 


background of processes which they are 


bringing into being. In addition, an administrative and 
financially minded class must be participants, the educa- 


tion of which is sufficiently broad toenable its members to 


recognize the sure foundations on which sound scientific 
conclusions rest. It is the object of the present publication 
to assist all three kinds of individuals... .” 

“ A study of the lists of contents of the twelve volumes 


which will eventually comprise the present series will 


indicate that an attempt is being made to acquaint the 
reader with the various stages of development of a manu- 
facturing process from the laboratory bench to the com- 
pleted factory. There are other matters, however, which 


must be taken into account if the factory is to be brought 
into successful operation and if its efficiency is. to be 
ensured. Consideration has been given, therefore, to 
Works Administration, 
The equally significant items of Works Safety and of 
Medical, First Aid and Welfare will 
included. Patent Legislation is another matter to which 
consideration has been given, an acquaintance with which 


may be all important to the « hemical engineer in the course 


Supervision and Maintenance. 


Services also be 


of exercising his wider functions.” 

The project is, therefore, an ambitious one, and one 
must admire the courage and enterprise of the editors 
in attempting it, and their success, so far, in obtaining 
contributions from well-informed authors. The completed 
work will be unique in chemical engineering literature, 
and will be found in all well-stocked academic and indus- 
trial libraries in this country, and many abroad. 

I suppose that any given well-informed reader whose 
somewhere in the chemical engineering 
novel information from 


Another third he will use as an aide- 


interests lie 


about one- 


spectrum will get 
third of the work. 
memoire in respect of matters with which he is familiar 
but not concerned frequently ; and the rest he will read 
as an expert—that is critically and suspi iously. It is 
essential, therefore, that the standard of all the contribu- 
tions should be high enough to compel the respect of the 
expert if he is to have confidence in the parts with which 
he is less familiar. 

A further requirement is that each contribution should 
be comprehensible to any potential reader (postulated by 
the editors to have reached degree standard in the pure 
or applied physical sciences). Finally, it is desirable to 
avoid too much overlapping, unevenness in presentation, 
diffuseness and the temptation to include more than 
can be adequately treated. 

No symposium has ever achieved perfection in these 
respects. A certain inconsequence and disjointedness, an 
invertebrate quality, is always apparent. For this reason 
one might speculate whether such talent and time as is 
available for writing technical books (all publishers will 
know how little there is) might not be more profitably 
used in producing a series of authoritative monographs. 
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In the present instance, might it have been more useful 
to promote a series similar to (but better than) the MeGraw- 
Hill Chemical Engineering Series ? However, a reviewer's 
task lies with the books he is given, and not with hypothe- 
tical alternatives 

The contributors to the present volumes, many of whose 
whok 


as far as I have checked 


names are well known, have on the done a good job 
Most of the information they give, 
it, is factually correct and the general level of presentation 
is good considering the difficulties involved in compressing 
the subject into the required compass. In a number of 
cases the treatment of topics is necessarily so condensed 
that it might have been better to take the working results 
for granted and refer to text-books for the theoretical 
bac kground The 
mics and of adsorption in “ The Preparation of Energy 
Balances” 


in the Solid State ” 


treatment of Second-Law thermodyna- 


, and of wave mechanics in “ Concepts of Matter 
are examples. 

On the other hand there are places where the editors 
to my mind, have allowed too much space to generalities 
For instance, 


with a rather low content of informatior 


the philosophy and methodology of pilot-plant work have 
been treated far too liberally, and the authors have plenty 
of space for stating the obvious, in sharp contrast to thos 
mentioned above 


of the 
is probably the most valuable part of these two 


The more technical sic section on pilot-plants 
volumes, 
although many of the topics dealt with under “ Design ” 
ilso be covered in Vol. IV 


and elsewher 


will presumably Design and 
Construction of Chemical Plant 


much of the thermodynamics dealt with in Vol 


Similarly, 


I is to be 


Volume 1 of Les 


Book reviews 


amount of solid 


tnnales du Genie Chimique has been received. 


covered again in Vol. IV. This tendency to duplication 
if each 


Nevertheless, 


is probably inevitable in a multi-volume work 
self-contained. 
it will add considerably to the total bulk of the published 


work, and it is to be hoped that the editors will do all they 


volume is to be reasonably 


can to co-ordinate the work of different authors. 


\ further danger which can be foreseen after looking 
at the first two volumes is that chapters may be included 
because expert authors are available rather than because 
their contributions are really relevant to the theme of the 
volume or of the work asa whole. For instance, E1iseNKLAM 
et al, 


Transpiration Cooling 


contribute an admirable chapter to Vol. II on 
this technique may one day 
find an application in chemical engineering, but the only 
applications at present are apparently to rocket motors 
and gas-turbine blades. A chapter on “ Porous Masses 


is followed by Industrial Uses of Porous 
Masses”. The 
the purification of coal gas, combustion in a 
and the 
ing as accounts of particular industrial processes, these 
the last) make little 


length in the 


* Important 
latter consists of three expert articles on 
fuel bed, 
filtration of water and sewage ; although interest- 
articles (except, to a certain degree 
reference to the prin iples discussed at 
earlier chapte r 

These 
the editors for the way in which they have largely over- 
difficulties of their task, and for the 
provided in their first 


criticisms are coupled with congratulations to 


come the inevitabl 


value they have 


two volumes. If they do as well or better with the remain- 


ing volumes, their work will become the vade-meceum 


of practising chemical engineers. , 


The series is intended to make 


available the results of work carried out in the Research Department of Institut du Génie 


Chimique of the University of Toulouse, 


Further volumes will be published annually 


Volume 1 (140 pp. 3.200 fr.) describes with detailed calculations the design of a process for 


the production of anhydrous sulphur trioxide. 
Directeur de Institut du Génie Chimique. 


The journal may be obtained from : 


It also contains a foreword by Professor Cathala, 


Editions Edouard Privat, 


14, Rue des 
Toulouse (France) 


Arts, 


vO 
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